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FOREWORD 

This report was prepared for the National Aeronautics and Space Ad- 
ministration Lewis Research Center under Contract NAS3-1 8544, 
Noise Addendum. Analysis of internal fluctuating pressure data (in- 
> eluding data from Phase I) and information about indirect combus- 

tion noise generated by interaction of temperature inhomogeneities 
with a full scale engine turbine pressure drop are presented here. 

« 

Appreciation is expressed for the contributions made by Mr. Ronald 
Huff of the NASA Lewis Research Center, and by Messrs. R. Roberts, 

I R. Dils, A. Peduzzi and G. Vitti of the Combustion Technology Group 

and G. F. Pickett of the Acoustics Technology Group at P&WA. 
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1. SUMMARY 


Internal measurements for the purpose of studying combustion noise characteristics were 
made on a variety of JT9D scale experimental combustion designs during the course of the 
NASA Experimental Qean Combustor Program, and partial results were reported in Refer- 
ence 1 . Tests were conducted on 90 degree combustor segments installed in a special bur- 
ner rig. Measurements were of two types: 1 . fluctuating pressure signals from transducers 
flush-mounted in the walls of the combustors, and 2. fluctuating temperature signals from 
an array of fast-response tl'.ermocouples located in the exit plane of the combustors. 

The objective of the pressure measurement program was to determine how the near field 
combustion chamber noise characteristics varied with combustor design and rig operating 
parameters, and to determine if certain operating parameters used to predict farfield core 
engine noise could also serve to correlate the internal measurements. It was found that a 
function of these parameters (Equation 9 of NASA TMX 71627, recommended by NASA as 
an interim prediction procedure for core engine farfield noise) provided a reasonably good 
description of changes in internal noise of a particular combustor as its operating parameters 
were varied. It was also found that the levels of internal noise of different combustor de- 
signs differed from each other by significant amounts at common operating conditions, and 
also that the average of the internal levels of the combustors evaluated at common operating 
conditions was appreciably higher than the corresponding predicted farfield value, as would 
be expected for information obtained from measurements in the near field of a noise source. 
Before these results were obtained it was found necessary to edit the internal pressure trans- 
ducer signatures to reduce the effects of local aerodynamic noise. 

The purpose of acauiring fluctuating temperature data downstream of the combustor was to 
obtain information relatlr" ' o what has been identified as “indirect combustion noise”. As 
distinct from direct combustion noise which is caused by local fluctuations in the heat re- 
lease rate, “indirect” noise is produced downstream of the burner zone when hot spots or 
temperature inhomogeneities in the flow encounter a sharp pressure change, such as is pro- 
duced by an engine turbine stage. No such sharp pressure drop existed in the operating of 
the combustor rig, so that no mechanism existed for the actual generation of indirect noise 
during the tests. The objective of taking the temperature measurements was to obtain infor- 
mation about spectra and length scales of the inhomogeneous flow which was then used as 
input to a computer program. Calculations were made, using performance characteristics of 
the JT9D high pressure turbine, to predict the indirect combustion noise that would be 
generated in a JT9D engine equipped with the experimental combustors. Details of the 
mathematical model used in the calculations are given in an Appendix to this report. 

It was found that the fluctuating temperature levels in one of the two combustors evaluated 
(vere appreciably higher than previously experienced with other, production engine combus- 
tors. Despite the high temperature fluctuations in this combustor calculated values of in- 
direct combustion noise were on the order of 10 dB lower than total core engine noise levels 
measured on a variety of engines and predicted by the previously referenced equation used 
in correlating internal pressure measurements. Therefore, it is expected that indirect com- 
bustion noise will not be significant in engines such as the JT9D operating with the experi- 
mental clean combustors tested in this program. 


2. INTRODUCTION 


As aircraft powerplant noise sources from the fan, jet, and turbine are reduced through 
improved engine design and use of sound-absorbing treatment, noise generated by the com- 
bustor could emerge as a significant limitation to further reductions. Compared to other 
powerplant noise sources, relatively little is known about combustion-generated noise. 

Noise measurements using combustor rigs have the advantage of allowing studies of the 
effects of varying operating parameters in a more independent manner than can be obtained 
in full scale engines. The NASA Experimental Clean Combustor Program involved emis- 
sions and performance tests on a variety of experimental combustors designed for reduced 
emission advanced CTOL aircraft engines. This program, conducted on a burner rig test 
stand, provided the opportunity to take internal measurements of the noise characteristics 
of several advanced combustor designs over a range of operating conditions. 

In Phase 1 of the program internal pressure measurements were taken on several combustors, 
and partial results were reported in Reference 1. During Phase II these data and further 
data on two other combustors were studied in greater detail. The objectives of this effort 
were to determine whether the internally measured noise varied with operating parameters 
in a manner similar to the predicted farfield noise, whether the internally measured noise 
levels were reasonably consistent with farfield noise . formation, and whether there were 
significant differences in the internal noise characteristics of different experimental com- 
bustors. 

A second type of investigation was also conducted during Phase II, which involved measure- 
ments of burner discharge temperature fluctuations. The purpose of this effort was to use 
these measurements to calculate and evaluate the significance of indirect combustion noise 
that would be produced by the interaction of the temperature inhomogeneities with a JT9D 
high pressure turbine. 
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3. COMBUSTOR EQUIPMENT ANC TEST FACILITIES 


During Phase I of the NASA Experimental Clean Combustor Program internal pressure data 
were taken during emissions tests on ten burner configurations which involved variations of 
three basic combustor types. Detailed descriptions are given in the ECCP Phase I Final Re- 
port, NASA-CR-1 34736, Ref. 3. In the Noise Addendum to Phase II both internal pressure 
and combustor exit temperature fluctuations were measured on two conflgurations, de- 
signated as H-6 and S-23, the features of which are described fully in the Phase II Final 
Report, NASA-CR-1 34969, Ref. 4. The principal design features of all the basic combustor 
types studied in the Noise Addendum are summarized below: 

Swirl-Can Combustor (N-Series) 

The swirl-can modular combustor design utilizes concepts developed by the NASA Lewis 
Research Center. This combustor incorporates a carburetor module array of 40 swirl cans 
in each of three circumferential rows, on a full annular basis. Each module is constructed 
with three major components: a carburetor, swirler, and flame stabilizer. The baseline con- 
figuration is shown in Figure la. Fuel is supplied to the modules through injection tubes 
centered in the carburetor-cans. The fuel-air mixture passes through the swirlers and is then 
burned. Dilution air enters the combustion chamber around hexagon shaped flame stabilizers. 
Fuel is staged to the outer row of carburetor cans at idle and to all the swirl cans at higher 
power. 

Six Swirl-Can Combustor configurations were evaluated during this program. These ’■e 
listed in lable 1 and described briefly in the sketches presented in Figures lb and Ic. lod- 
ificatiens were confined to diffuser inserts, in an attempt to provide more uniform difiuser 
flow, and carburetor module changes to alter flameholding characteristics and carburetor 
fuel-air ratio. In one modification, configuration N-1 1 , inner and outer liner dilution holes 
were provided. 

Staged Premix Combustor ( P-Series) 

The design of the Staged Premix Combustor is based on previous single and two-stage pre- 
mix combustor developments conducted at Pratt & Whitney Aircraft. Figure 2a presents the 
baseline design involving use of two comoustion zones, denoted as pilot and main. Each 
zone has a separate fuel injector array, premix passage, flameholder, and combustion volume, 
as may be seen in Figure 2a. .Two zones are employed with the intent to facilitate emission 
control over the range of engine operation more readily than can be accomplished with a 
single zone combustor. The pilot zone is sized to accommodate all heat release for idle 
operation and is designed for optimum emission control at idle. With appropriate fuel man- 
agement the combination of both pilot and main systems can be optimized for the high 
power operating regime. 

Two configurations of the basic Staged Premix Combustor were tested in Phase 1 as indi- 
cated in Table 1. Both incorporated slotted pilot and main flameholders, as shown in 
Figure 2b with no liner dilution air. This produced a very lean bulk premix passage equiva- 
lence ratio. Configuration P-8 additionally incorporated 2X main zone fuel source density. 
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TABLE I 


COMBUSTOR DESIGN MODIFICATIONS EVALUATED 
IN ACOUSTICAL PROGRAM 


Combustor 

Scheme 

Configuration 

Figure Ref. 

SwirhCan Combustor 

N-7 

Diffuser Screen; 1.27 cm recessed 
swirlers 

lb, Ic 


N-8 

"V” gutter trip; non-recess:, swirlers 

lb, Ic 


N-9 

Diffuser screen; outer swirler flame- 
holder design. 

lb. Ic 


N-10 

Diffuser screen; reduced swiri can air. 

lb. Ic 


N-11 

Diffuser screen; ID and OD dilution 
air. 

lb. Ic 


N-12 

Diffuser screen; outer swirier fiame- 
holder design; pressure atomizing 
nozzles. 

Ib, Ic 

Staged Premix Combustor 

P-7 

No dilution air; slotted flameholders 
pilot and main. 

2a. 2b 


P-8 

No dilution air; slotted flameholders 
pilot and main; 2X main fuel nozzles. 

2a. 2b 

Swirl Vorbix Combustor 

S-8 

No dilution air; high flow main 
swirlers; pressure atomized pilot fuel 
nozzles. 

3 


S-9 

15 percent pilot swirlers blockage; main 
fuel nozzles aimed downstream. 

.7 

Swirl Vorbix Combustor 
(Phase II) 

S-23 

Pressure atomized pilot nozzles with JT9D-70 
swirlers; reduced main zone swirler airflow; 
increased pilot volume; dilutic i air in both 
pilot and main zones. 

4.5 

Hybrid Combustor (Phase II) 

H-6 

Pressure atomized hollow cone primary 
nozzles; main zone swirler airflow increased 
40%; no dilution air. 

6, 7 
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Swirl Vorbix Combustor (Series) 


The Swirl Vorbix (Vortex Burning and Mixing) Combustor also employs two combustion 
subsystems, as shown in Figure 3. The pilot zone is a conventional swirl stabilized, direct 
fuel injection design and is sized to provide all heat release at high efficiency during idle opera- 
tion. At higher power conditions, main zone fuel is introduced through fuel injectors located 
downstream of the pilot combustion zone. Main combustion and dilution air is introduced 
into the main burning zone through 60 swirlers located on each side of the combustor (on 
an annular basis). The hot exhaust product - fuel vapor mixture is entrained by the swirling 
air jets and partial ptemixing occurs prior to auto ignition of the main zone fuel. The main 
combustion process, in which most of the fuel is consumed at high power conditions, pro- 
ceeds rapidly at the interface between these jets and the surrounding fuel rich pilot jet mixture. 

Two configurations of the Swirl Vorbix Combustor were evaluated in Phase I, S-8 and S-9. 
Configuration S-8 differed from the baseline design in that pressure atomizing fuel nozzles 
were substituted for the aerating pilot zone nozzles, and large swirlers were installed at the 
main throat location. Th.-se had somewhat more than twice the effective flow area of the 
baseline swirlers, and necessitated the elimination of all aft liner dilution holes in order to 
maintain pressure drop. S-9 incorporated pilot swirler blockage rings to increase pilot zone 
equivalence ratio, and modified the orientation of the main fuel injectors. Tlie downstream 
injection of main fuel was intended to minimize combustion prior to mixing with the main 
swirler airflow. 

Swirl VorbLx Combustor S-23 (Phase III 

The basic Swirl Vorbix Combustor design used in the Phase 11 program (S-1 1) differed from 
the configuration evolved during Phase I in that the length of the pilot burner was increased 
by 3.8 cm (1.5 inch), while the main burner length was reduced by this amount. These 
changes were made to improve emissions characteristics. The specific model on which ac- 
oustic tests were conducted, S-23, differed from the basic Phase 11 Swirl Vorbix design in 
using JT9D-70 pilot swirlers and having reduced main burner swirler flow. The pilot zone 
volume was increased about 30 percent, and dilution air was added to both pilot and main 
zones. The basic Phase 11 design features are illustrated in Figures 4a, b, and c. Additional 
details applicable to configuration S-23 are showm in Figure 5. and a more extensive descrip- 
tion is found in the companion main ECCP Finai Report NASA CR-1 34979, Ref. 4. 

Hybrid Combustor H-6 (Phase II) 

TTie specific Hybrid Combustor, H-6, on which noise data were acquired in Phase II was a 
modification of the basic Hybrid Combustor, H-1 , shown in Figures 6a and 6b. The com- 
bustor employs two burning zones that are separated radially and axially. The design basis 
for each zone evolved from two potentially good concepts for high and low power emission 
control that were developed in Phase I. The pilot zone duplicates the geometry of the Premix 
Combustor P-3 pilot which exceeded the low power emission goals in Phase 1. The main burn- 
er zone design was based on the Swirl Combustor N-9, which demonstrated significant reduc- 
tions in high power emissions reduction. 
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Fuel is introduced into the pilot premix passage through thirty pressure-atomizing fuel inject- 
ors (on a full annular basis), where it mixes with air and is partially vaporized. This mixture 
then passes through holes in the flameholder into the pilot combustion zone. The main burn- 
er zone has an amngement of concentric, counter-rotating swirlers, located downstream and 
displaced radially from the pilot zone. Fuel enters through low pressure drop iryectors in the 
center of the swirlers. 

The specific Hybrid Combustor, H-6, evaluated for noise characteristics in Phase II incorpor- 
ated modifications that increased the outer swirler airflow about 40 percent and blocked the 
inner liner dilution air. Details of this configuration are shown in Figure 7. 

Test Facility and Combustor Rigs 

The cv mbustor noise tests were conducted on X903 stand, a high pressure combustion com- 
ponent test facility located at the Pratt & Whitney Aircraft Experimental Test Department, 
Middletown, Conn. This facility was fully equipped with the necessary inlet ducting, exhaust 
ducting, controls and instrumentation for conducting component tests at simulated engine 
operating conditions. For testing, the combustor test rig containing the combustor was instal- 
led at the inlet and outlet flanges of the test stand ducting. Heated non-vitiated compressed 
air was supplied to the test rig and the discharge pressure was controlled by valves in the ex- 
haust ducting, allowing operation over a range of test conditions. 

Different complete 90 degree sector combustor rigs were fabricated for use in this program. 

A JT9D-7 engine diffuser/burner case was sectioned to provide the same strut oriei’tation in 
all rigs. Each case was modified to provide the fuel support pads and instrumentation bosses 
required by the particular combustor design. Stand duct work and flanges were also fabricated 
to mate the rig to the facility. A schematic diagram of the test rig, showing the location of 
the performance instrumentation is given in Figure 8, together with a diagram of the overall 
test facility. 

Instrumentation was provided to measure the variables required to determine combustor per- 
formance, such as air flow, fuel flow, fuel temperature, static pressure and total temperature 
at the combustor inlet, and total temperature and pressure of the gases at the exhaust of the 
combustor. Results of these measurements were used in correlating noise with operating 
parameters. 



4 . ACOUSTICS INSTRUMENTATION AND TEST PROCEDURE 


Pressure Measurements 

Combustor internal pressure measurements were made using Kistler Model 603AI dynamic 
pressure transducers. The transducers were installed in Kistler water cooling jackets, model 
616M13S, in order to protect the transducer from the extreme temperatures. The transducer 
installation is shown schematically in Figure 9. Additionally, in Phase I an inlet transducer 
was located on the ID wall of the combustor inlet duct, 0.299m upstream of the diffuser case 
forward mounting flange. 

The sidewall combustion compartment transducers were located in the same relative locations 
in all combustor configurations. The locations were in the combustor sector sideplates, cen- 
tered radially approximately 0.08S9m forward of the combustor rear mount flange. Details 
of the combustor sidewall transducer installation are shown in Figure 10. Due to procedures 
required for installing the combustion compartment transducers, it was necessary to install 
them as the combustor was being assembled. This required at least two separate sets of 
transducers allowing one combustor to be built up while another was at test. 

Ifee of the pressure transducers requires determination of two different characteristics: ( 1 ) 
sensitivity and (2) frequency response. “Sensitivity” is the output voltage per unit sinusoidal 
input pressure at a fixed reference frequency. “Frequency response” is the relative change in 
output voltage as the frequency of a constant amplitude input pressure is varied. During 
Wiase 1 it was determined (Ref. 1) that the frequency response was essentially constant over 
the range 30 to 5000 Hz. Small deviations were noted for some of the transducers and were 
applied in the data reduction procedure, in Phase II for tests of Configurations H6 and S23 
it was decided that these corrections were sufficiently small to ignore. 

The transducer sensitivity was determined by use of a system involving a pressure source pis- 
lonphone, the calibration of which is traceable to NBS standards (Ref. 1), a previously-calibra- 
ted reference transducer, and the subject transducer. Both transducers employ charge ampli- 
fiers in the system. The reference transducer and charge amplifier system (NBS traceable) is 
subjected to the pistonphone output, operating at a sound pressure level of 170 dB and a 
frequency of 1 5 Hz, and the system output noted. The reference transducer-charge amplifier 
system is then removed from the pistonphone assembly and replaced by the particular Kistler- 
amplifier being calibrated. By setting and locking the subject charge amplifier gain to produce 
the output voltage recorded for the reference system, sensitivity of the subject transducer- 
amplifier is set equal to standard. No further correction procedures are required in data acqui- 
sition of data reduction. After completion of testing, the Kistler systems were re-checked 
and were found to be operating properly. 

The data recording system consisted of: 1 ) the pressure transducers and their power supplies, 
and 2) the signal conditioning and recording equipment, which provided calibration and 
monitoring instrumentation, switching capability, and analog magnetic tape recording capa- 
bility. Recordings were made at a tape speed of 30 inches per second, using the FM mode 
of a Precision Instrument, Model 200 Recorder. Figure 1 1 is a schematic diagram of the data 
recording system. Recording system response calibration signals were provided by a local os- 
cillator and random noise generator. Gain settings established on each channel at the pre- 
amplifier, to obtain optimum signal to noise recordings, were noted on the recorder log 
sheets and announced on the tape voice channel. 
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The data analysis system is shown in Figure 12. The recorded pressure signals were fret^ ^ency 
analyzed using a set of 24, one-third octave band filters having geometric mean frequencies 
extending from 50 Hz to 1 OK Hz. These filters complied with the filter characteristics recom- 
mended in lEC Publication 225. The system was set up to provide a readout of one-third 
octave levels for each transducer, time averaged over a minimum time period of 30 seconds 
for each operating condition. Following this analysis, the data were stored on digital mag- 
netic tape for input to the computer. Recording response and transducer frequency response 
corrections were applied to the data in the computer from SO Hz to SK Hz. 

Selected narrow band analyses were performed o*' a high speed spectral analysis system pro- 
viding wide-range, high-speed digitization of the narrow band analysis analog data. Resultant 
spectra were stored and averaged in a local memory before being plotted as narrow band 
spectrum plots. 

Temperature Measurements 

Through use of fine wire elements, carefully butt-welded at the junction, thermocouples have 
been developed that operate reliably in engine combustor environments, and with electronic 
compensation, will respond to dynamic temperature fluctuations in the stream in the fre- 
quency range from 0 to 1000 Hz. 

In this program, an array of 16 such elements was used to measure the dynamic temperature 
field of the discharge from two combustor designs. Details of the thermocouple elements are 
described here together with features of the 16-element array installation, and the recording, 
playback, and signal processing techniques used in this investigation. 

These Type S thermocouple elements were formed from pairs of .00762 cm (.003 inch) dia- 
meter wires, one of platinum and the other of platinum with 10% rhodium. The junction 
between each wire pair was formed by laser beam butt-welding. The resulting wire element, 
with the weld at midpoint had a span of .23 cm (.09 inch), and was spot welded at each end 
to .025 cm (.010 inch) diameter support posts. These posts are .5 1 cm (.2 inch) long stripped 
ends of standard Type S thermocouple wires sheathed with MgO insulation. 

Four such thermocouple elements, spaced 1 .08 cm (.425 inches) apart were installed in a 
rake assembly which was installed in the combustor rig to measure temperatures at four 
radial locations in the annular burner discharge station. Figure 13 shows details of the ther- 
mocouple elements and the rake construction, which incorporated a steam cooling jacket for 
increaang the durability of the thermocouple wire insulation and the rake structure itself. 
(The thermocouple element is not affected.) The fine wire elements can survive environmen- 
tal temperatures on the order of 2500 to 3000° F long enough for most test purposes, but are 
subject to the hazards of overtemperature transients during ignition and other changes in 
setting operating points. To guard against such premature failures, the rake assembly was 
fitted with another jacket system to provide jets of screening air for the thermocouple ele- 
ments during the setting of rig operating conditions. These jets were arranged to not actually 
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impact the wiie element but rather to act as aerodynamic windscreens by diverting the bur* 
ner efflux locally away from the element. After stabilization of the operating point condi- 
tions the screening air supply was turned off prior to talcing data. Four such r^es, each 
mounting four fine wire thermocouple elements were mounted in the discharge plane of the 
annular combi stion chamber segment. This array is shown in Figure 14. The circumferential 
separatr >n between elements on adjacent rakes varied between about 7 to 8 times the 1.08 cm 
radial s^par. <' on. 

Details of th(; heat transfer characteristics that determine the dynamic response function of 
these thermocouples may be found in References 9 and 10. Very briefly, the wire element 
acts as a low frequency bandpass filter for harmonic components of the unsteady input tem- 
perature. The “cutoff frequency” depends upon local steadystate heat transfer characteris- 
tics of the wir ^/gas system and in practice has to be determined (usually in terms of the cor- 
responding me constant”) for each operating point. Roughly, the basic wire cutoff freq- 
uency is very ow - on the order of 10 Hz. Above this frequency the output voltage response 
drops at a rate approaching 6 dB per octave. To obtain information on temperature fluctua- 
tions boyond about 10 Hz, electrical compensation must be provided. From previous 
experience a frequency range of zero to 1000 Hz was selected as providing a balance between 
true teniperati-re fluctuation information and unwanted amplification of electrical system 
noise. Compensating amplifiers have been developed with the properties of balancing the 6 
dB per octave drop above thij bare wire cutoff frequency with a transfer function such that 
the overa!' system frequency range is extended by a factor of about 100. At high frequencies 
the overai; system response (including the compensating amplifiers) falls at a rate approaching 
12 dB per octave. In use, the network has an adjustment which must be set to the wire time 
constant i.i order that compensation :larts at the correct frequency. This time constant is 
determine^ by making a simple calculation involving knowledge of the steady state gas pro- 
perties at t ie operating point. Figure 1 5-a is a diagramatic representation of the thermo- 
couple and compensating am*^ !ifier system. Figure 1 5-b gives the voltage frequency response 
characteristic of the overall system, and also shows a sample uncompensated thermocouple 
response. Tire system “cutoff’ frequency, f„, (defined as that frequency for which the com- 
bined system response is down 6 dB) in the response function expression has been set from 
previous experience at 1000 Hz. (While it might appear desirable to modify the network to 
provide still higlrer frequency range, such an alteration would have the effect of amplifying 
unwanted electronic system noise; the cutoff frequency of 100 Hz has been fixed on the 
basis of considerabie development experience with these instrumentation systems.) 

In tests employing relatively few thermocouples, each thermocouple is provided with a com- 
pensating amplifier in the ape recording data acquisition process. However, because this 
test program used a large (16) number of thermocouples, an adequate number of compen- 
sating networks cc..ld not be provided for recording. Instead, the uncompensated signals 
were amplified ....J recorded directly, using the Precision Instrument FM tape recorder pre- 
viously described in connection with the pressure transducer signals. (The 40 dB dynamic 
range of f tape recorder allows valid recording to be made without use of compensating 
amplif us.) Two records were taken at each operating point, selecting the pressure trans- 
ducers and some of the thei mocoupies for one record and the remaining inputs for the 
second recording. 
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During playback of the thermocouple information, either one or two compensating ampli- 
fiers were inserted between the tape recorder and the signal analyzer. Two essentially dif- 
ferent types of signal processing were used for the thermocouple information. One-third 
octave band voltage spectra were obtained using a single compensating amplifier and the pre- 
viously described wave analyzer. Values of overall root mean square voltage fluctuation were 
obtained from a B&K 2425 true rms voltmeter. To obtain the corre^onding rms temperature 
fluctuation the ASTM Star, iard Temperature-Electromotive Force Tables for Thermocouples, 
E-230-72 were employed. 

The table for type S thermocouples was entered at the mean steady state discharge tempera- 
ture of the operating point. The local sensitivity computed from the table at this point, 
when multiplied by the rms voltage fluctuations, gave the corresponding rms temperature 
fluctuation. This sensitivity was essentially constant over the test operating range with an 
average value of about 0. 1 5 degrees Fahrenheit per microvolt. 

The other mode of signal processing involved correlation, using a Saicor Model 43A correla- 
tor. For autocorrelation of selected thermocouple signals a single compensating amplifier 
was inserted between the tape recorder and the correlator. For cross correlation analysis of 
selected pairs of thermocouple elements two compensators were provided. The purpose of 
the correlation processing was to determine axial and transverse length scales of the tempera- 
ture inhomogeneities in the combustor discharge gas stream. Details of how these scales were 
obtained from the correlograms are presented in the Results section of this report. 

Test Procedures 

In Phase I of the program dynamic pressure measurements were taken for the combustors 
operating at simulated idle and takeoff conditions. Corresponding cold flow (fuel ofQ condi- 
tions were also recorded. In Phase 11, tests on Configuration H-6 also included simulated ap- 
proach operation and the corresponding cold flow point. Tests on Configuration S-23 further 
added operation at climb conditions. Temperature data were taken only during the H-6 and 
S-23 tests. 

Combustor performance data were also taken to provide information for correlating the fluc- 
tuating pressure and temperature information with steady-state performance parameters. 
Measurements were made to obtain airflow, fuel flow (primary and secondary fuel systems), 
combustor inlet and total static pressure and total temperature, inlet mass flow parameters, 
combustor liner pressure drop, and ideal average combustor exit temperature. 

In Appendix A-1 are tabulated those parameters of greatest significance from the standpoint 
of identifying the run conditions and providing a particular correlation between noise and 
combustor performances. These tables arc supplemented by acoustic power level listings 
that are defined in the Results section of this report. They are listed in Appendix A-1 to 
provide a convenient reference location. All points recorded in Phase II are given in these 
tables. However, detailed analysis of the data acquired in Phase 1 was limited here to certain 
points for cc ifigurations N-8, N-1 1, P-8 and S-8, and consequently only those corresponding 
points are tabulated. Performance parameters for the other tests conducted in Phase 1 are 
given in Reference 1. Further details of the construction, performance, and emission features 
of the Phase II configurations may be found in the FCCP Phase H Final Report, NASA-CR- 
13496‘). 
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5. RESULTS AND DISCUSSION 


Because of the significant differences in basic phenomena, measurement systems, and analyti- 
cal procedures, between the internal pressure and combustor exit temperature investigations 
conducted in this program, the presentation and discussion of the results is divided into two 
sections: first, the pressure data, and last, the temperature information. Furthermore, back- 
ground material necessary to understanding the specific quantified results is presented in 
subsections entitled, “Method of Analysis” which precede the specific “Results” subsections. 

5.1 PRESSURE INFORMATION 

For convenience in discussing the pressure information, the objectives given in the Introduc- 
tion to this report are repeated here: These objectives were to determine the characteristics 
of internally measured noise for a variety of experimental combustors and to note whether 
there were significant differences in the overall noise among the configurations. It was also 
sought to study how variations in operating parameters affected the noise, and in particular 
to determine whether this relation could be quantified by a simple expression recommended 
for use in predicting farfield core engine noise. 

5.1.1 Method of Analysis 

Corrected Sound Pressure Levels 

For each set of runs at various degrees of burr jr temperature rise obtained at the same mass 
flow and burner entrance pressure and temperature, a run was made with no burning and is 
referred to as the corresponding “cold” point. Cold points were obtained for simulated idle, 
approach and takeoff conditions for configuration H6, and in addition, at climb for con- 
figuration S23. By examining a set of spectra for a particular condition, e.g., idle, together 
with the cold point in that set (Appendix A-3) it is evident that the effect of burning on in- 
creasing the noise is generally limited to a frequency range from about 1 00 to about 4000 
Hz, and is usually evident over a narrower range. A major portion of the pressure transducer 
signal is due to noise of aerodynamic origin rather than unsteady combustion and does not 
constitute a valid measure of the process under study. To obtain a better measure of the 
noise associated with burning, therefore, requires that the “as measured” spectra be cor- 
rected in some manner to de-emphasize the aerodynamic noise contribution. The method 
selected here involved subtracting (loprithmically) the corresponding cold flow spectrum 
from each of the spectra obtained with the burner operating. Further, where necessary, no 
contributions to burner noise were assigned below 1 25 Hz or above 4000 Hz. By summing 
(logarithmically) these “corrected” 1/3 octave band levels in the 125-4000 Hz range, a 
“corrected” overall sound pressure level, OASPL^.^,, was obtained. In this manner a single 
number rather than a spectral distribution function could be used as a measure of the com- 
bustion noise signal from the transducer and employed in subsequent calculations. The 
uncorrected overall value, OASPL, was also employed in a limited way, as will be described 
presently. Sample spectra, illustrating the correction method are given in Appendix A4 for 
H6 and S23 and 4 configurations selected from Phase I for analysis. 
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Power Level Calculation. PWL and PWL 


Corresponding to values of OASPL and OASPL^.^,, the acoustic power levels, PWL and 
PWL^„, were obtained from the logarithmic form of the expression for power: 

P = A ( 1 + M), where 
pc 

P = acou tic power 

p^ = mean square pressure 

p = fluid density 

c = speed of sound 

A = section area 

M = Mach number, 

where all quantities are evaluated in the plane of the transducer, are in mks system units. In 
the calculations A was taken as 4 times the 90 degree rig sector area to obtain representative 
values for a full-scale complete annulus. The Mach number in the plane of the transducers 
was always less than 0.09 and was therefore ignored since there are much larger uiicertainties 
in the acoustic power determination as explained twlow. 

In the mks system a reference PWL of Iff* ^ watts is used. To obtain P’.VL values in a system 
where Iff*^ watts is taken as reference (see following section on PWLp^^j ) 10 dB must be 
added. The following expressions give the acoustic power level (re 10-13 watts) fyjj 

scale engine where Mach number is ignored. A is the 90“ sector rig area (the correction to 
the full annulus area is included in the constant), and all units (except PWL as noted above) 
are mks units; 

PWL = 42 -HO logA - 10 log pc -i- OASPL 
•^1-cor = above -i-OASPL^„, 

Near I 'ield, Transverse Modes, and Transmission Effects 

The preceeding equation for acoustic power is correct only when the acoustic particle velocity 
and pressure are in phase, and when the direction of wave propagation is locally normal to 
the e'l ment of surface, A. Both restrictions are necessarily violated in the very near field of 
a source distribution such as was the case with the transducer measurements in this program. 

For a discrete frequency or band-limited coherent wave structure, if the particle velocity, u, 
could be determined, the quantity ^/pc should be replaced by p cos 4> where 




0 is the phase angle between pressure and particle velocity. Therefore, by employing p^/pc 
to compute “power” the effect of phase angle, cos 0, or “power factor” is necessarily sup- 
pressed. For this reason it can be seen that the values of “power” as computed here must be 
higher than would be measured at some distance from the source region, where the acoustic 
pressure and particle velocity are actually in phase. Due to the construction of the burner 
test rig it was impossible to position microphones downstream of the combustor in a region 
where they would not be impacted by disturbed flow or subjected to other undesirable side 
effects. 

These limitations should not be construed as an argument against the use of combustion 
chamber pressure transducers, for in an engine there is no other location available for internal 
measurements that is free from equivalent complications. It would be highly desirable, how- 
ever, in future combustion noise experiments to provide the test apparatus with a downstream 
duct location sufficiently far from the combustor so that near field effects could be avoided 
and such that flow, valve, and vibration noise effects were minimal. This location would be 
used to obtain true acoustic power measurements which in turn could be used to correlate 
with the combustion chamber readings. 

Another limitation of the data (which would also apply to a downstream duct location free 
from the above mentioned near field effects) is that one or two transducers in general cannot 
adequately define the acoustic field structure in a duct whether they are close to or well 
separated from the source region. Under conditions where Trequency is low enough such that 
only the plane wave can propagate, a single transducer can define the field (if there are no 
end-reflections). However, where the range of frequencies is such that higher transverse duct 
modes are above cutoff, the acoustic pressure distribution in a transverse plane of the duct 
will not be uniform (as it is for the plane wave only). In such cases a plurality of transducers 
must be employed in an array and their outputs must be specially processed in order that the 
mode structure and power flux can be computed properly. Tlie dimensions of the test com- 
bustors used in this program and the general temperature range of operation were such that 
the first transverse mode could propagate for frequencies greater than about 500 Hz. There- 
fore, in addition to the near field effects described previously, the presence of transverse duct 
modes makes the estimation of acoustic power by means one or two transducers subject to 
errors of uncertain magnitude. 

The expression P = p^ A/pc, used here to calculate power, is known to overestimate the 
true farfield power due to both near field effects and the presence of transverse modes. The 
near field “power factor” effect has been explained. It can also be shown that a transverse 
mode (in a rectangular duct) may be synthesized by a pair of plane wave trains propagating 
at an angle, /?, with respect to the duct axis. If the duct cross section area, A is used, the 
power flux through this^ea will involve the component of the inclined wave intensity along 
the duct axis or = (p^ A/pc) cos where the m-subscript refers to a particular trans- 
verse mode. Since the calculation of “power” in this program has had to be based on assuming 
both power factor and direction cosines to be unity, the values thus obtained will be over- 
estimates of the correct figures. In an engine, due to the effects of passage through the turbine 
stages, stream turbulence, and terminal nozzle impedance, not all of the “true” combustion 
acoustic power is radiated to the far field. These internal losses are a tliird cause for the 
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“power” as obtained from the internal transducer readings to overestimate the combustion 
acoustic power that will, in fact, radiate to the far field. 

Aerodynamic Noise Check Procedure 

In the above section dealing with corrected sound pressure levels it was stated that the 
transducer signal spectra indicated significant contributions from aerodynamic noise. (When 
the burner was ignited, only certain third-octave band levels increased.) For verification steps 
were taken to determine if: (1) the overall power level, PWL, could vary with rig operating 
conditions in a way that was consistent with an aerodynamic noise source; and (2) that the 
corrected power level, PWL^^^, had in fact been freed from this influence. 

Noise of aerodynamic origin may be characterized in the far field by a super-position of 
source distributions that are referred to as monopole, dipole, and quadrupole in behavior. 

If one of these source distributions dominates, the far field acoustic power will vary as the 
stream velocity raised to an exponent that is 4 for a monopole distribution, 6 for dipoles, and 
8 for a quadrupole source field. However, regardless of the far field behavior, the pressure 
fluctuations in the near field of any aerodynamically generated noise will vary as pV^ where 
V is a characteristic mean velocity. 

The “power” P, calculated from the square of the near field pressure, is thus proportional to 
. Therefore, in logarithmic units, this overall power level, PWL, will equal 4(10 log v) 

+ 20 log p + a reference constant if the source is aerodynamic. 

When the quantity (PWL - 20 log p) is plotted against 10 log V a line with a slope of 4 results. 
Hence if a plot of these variables is made for a set of test points, and a line with slope of 4 
results, there is strong evidence that the pressure transducer is sensing the near field of an 
aerodynamic source distribution, rather than (as is desired in this program) noise generated 
by unsteady combustion. Accordingly, to test whether PWL was possibly dominated by 
flow noise, and esi^ecially to confirm that the procedure for computing the “corrected” 
quantity, PWL^.^^^ had removed this influence, plots of the type described were made for tests 
conducted on Configurations H6 and S23 using average values of stream density and velocity 
in the discharge plane of the combustor, station 5. This information is presented at the be- 
ginning of the “Results” subsection below. 

Correlation of Pressure With Rig Operating Parameters 

For several reasons, listed below, the task of correlating the internally-measured pressure 
fluctuations with rig operating parameters was performed by employing equation (7) of 
Reference 2, “NASA TMX-71627. Interim Prediction Method for Low Frequency Core Fn- 
gine Noise”. When the quantities in the referenced equation are expressed in the notation 
used here, Fq. (9) becomes: 
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= 56.5+ 10 log [(Tj -T 4 ) 

P T, 

o 4 

where PWLpjgj = overall far-field acoustic power level, dB re 10*’^ watts 
= engine airflow (Ibm/s) 

T * Total temperature (°R) 

P = Total pressure (Ib/ft^) 

Station 0 = atmospheric 

4 = combustor inlet 

5 = combustor exit 

In the rig tests, a 90 degree sector of the annular combustor was used, so to obtain corres- 
ponding engine predictions 6 dB should be added if rig airflow is used, making the constant 
62.5 instead of 56.5. 

Tlie method used here for correlating measurements with operating parameters consisted of 
computing PWLp|.gj (called predicted power level) for the operating points and comparing 
the measured and predicted power levels. As previously explained, measured power employed 
a correction for cold flow noise, and the resulting quantity is called PWL^^^. 


The expression for PWLp^^.^ as a correlator function for the internal pressure measurements 
was selected for these reasons: 


1. Eq. 9 is reported (Ref. 2) to do a good job of predicting far field low frequency 
core engine noise. Tlie constant differs by 16 dB depending on engine type (e.g., 
turbojet, turbofan, turboprop) but not with engine size. 

2 . Item 1 is an important reason for determining whether this equation (with a dif- 
ferent constant) also is a good predictor of combustion noise measured internally. 

3. Eq. 9 has the virtue of being simple to use. 

For the above reasons, the combination of rig operating parameters given in Reference Eq. 9 
was used as the basis for correlating noise. The mechanics of performing this correlation 
was extremely simple: at each operating point, PWl.p,gj was calculated by ihis equation and 
used as the abscissa of a point having ordinate given by (measured) PWL^.^,. If the resulting 
plot of PWL^^j vs PWL scatters about a line with slope of unity, PwL^j is a good pre- 
dictor or “correlator” for the measured power. To evaluate the degree of “goodness” the 
standard error of estimate (standard deviation using N-2 degrees of freedom) was obtained 
and compared with the corresponding deviation of the data points from a horizontal line 
through the centroid. The latter quantity is a measure of the null hypothesis - - that the 
measured quantity PWL^^ does not depend on PWL 
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5.1.2 Re$ults 


The following results were obtained by applying the above procedures to the data: 

Pressure Data Presentations 

Pressure data tabulations and plots for tests conducted on configurations H6 and S23 are 
presented in Appendix A, and corresponding information for the configurations tested in 
Phase I is contained in Reference 1, Appendix A-1 tabulates the ECCP rig aerodynamic per- 
formance data and may be used to identify the operating conditions corresponding to the 
point numbers for which acoustic data were taken. In Appendix A2 are found tabulations of 
1/3-octave band and overall sound pressure levels for both left and right hand side transducers. 
Machine-generated plots of 1/3 octave band spectra for the right hand side transducer are 
contained in Appendix A-3. No plots are given for the left-hand transducer because this in- 
formation was not used in the subsequent analysis. Although there was not evidence of 
transducer malfunction during the recording of data for tests on H6 and S23, nor during post- 
test calibration, examination of the tabulated and plotted spectral information showed serious 
discrepancies at several operating conditions between the two microphone spectra. Order of 
magnitude studies were made of possible extraneous transducer response to vibration, heat 
fluctuations, aerodynamic noise and local flow impact noise. Based on these studies it ap- 
peared that the left microphone frequently indicated levels that were inconsistent with changes 
in rig operating conditions, and was probably more affected by local aerodynamic noise effects 
than the right hand transducer. 

Aerodynamic Noise Checks 

Estimation of possible microphone signal domination by aerodynamic noise sources rather 
than combustion noise was performed for configurations H6 and S23 by computing overall 
power level corrected for density, obtained by subtracting 20 log pj from PWL, and plotting 
the result versus 10 log Vj. Figures 16 and 17 present this information. Least squares fits 
give lines with slopes of 4. 1 for H6 and 0.7 for S23, as shown on the figures. The slope for 
H6 is uncomfortably close to 4, indicating probablp response to noise of aerodynamic origin. 
Spectra for H6 compared with cold flow spectra show that the high frequency portions 
(above 4K Hz) are generally similar and of large amplitude, confirming that the overall sound 
pressure level (from which PWL is obtained) is a poor measure of noise due to unsteady burn- 
ing. In spectra for configuration S23 the influence of burning predominates over a wide freq- 
uency range, accounting for lack of V** dependence. 

Therefore, overall power level as such cannot be used with confidence as a measure of com- 
bustion noise, and a correction procedure must be applied. The procedure used here to ob- 
tain a corrected power ’evel, PWL^.^,,, involved logarithmically subtracting the cold flow noise 
spectrum from the spectra obtained with burning at common values of combustor mass flow 
and entrance pr assure and temperature and truncating the result at 4K Hz. When this cor- 
rected value, PWL|,jjj, is used instead of PWL the results are shown in Figures 18 and 19. 

Slopes of the least squares lines of 5.7 and 2.9 for H6 and S23 suggest that aerodynamic noise 
has been eliminated as a dominant influence. 
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Alternative procedures exist for emphasizing the combustion noise component of the internal 
transducer signal. In full scale engines rather than rig tests no cold flow data acquisition is 
possible. There is sufficient general evidence in the combustion noise field that noise due to 
unsteady burning is limited to relatively low frequencies. Therefore if the spectra are arbitrarily 
truncated at an upper frequency limit such as the 4K Hz employed here the combustion 
noise contribution will be emphasized. 

Configuration H6 Results 

As described in the subsection on method of analysis, study of the relation of combustion 
noise to rig operating parameters was confined to a particular function of the parameters de- 
fined by Eq. 9 of TMX-71627. Values of acoustic power predicted by this equation, called 
PWL were computed for each operating condition and used as the x-coordinate of a 
point, the corresponding y-coordinate being obtained from pressure transducer mea- 

surement. 

If the points thus obtained define a line with slope of unity, Eq. 9 is established as a perfect 
predictor of internally measured acoustic power. The standard deviation or standard error 
of measurement of the actual set of data points is a measure of systematic or random de- 
viations from this predictor, large values indicating poor correlations. To provide a reference 
value against which this deviation can be compared the standard deviations from a horizontal 
line through the centroid of the set of points may be used. This procedure is comparable to 
testing the “null hypothesis” that measured noise is completely uncorrelated with the com- 
bination of operating parameters given by PWEp^^j. 

Figure 20 shows values of PWL^,^,j versus PWLp,^j for the 15 operating conditions tested. 

Lines with unity and zero slopes are drawn thru the centroid of the distribution. Standard 
error of estimates for the data with respect to Eq. 9 predictor is 2.75 dB (and 4 dB for the 
flat line), indicating that Eq. 9 is a reasonably good predictor, 95% of the data falling within 
±5.5 dB of the predicted values. It is worth noting that the constant difference, about 16 dB, 
between the slope 1 line through the data and the line through (0, 0) is compatible with near 
field and far field differences and to reported differences in the constant among far field data 
for different engine types (Reference 2). 

limited tests in Phase I indicated that once the burner was ignited at an operating pressure, 
temperature and flow, changes in temperature rise seemed to have small effects on the spectra. 
iTiis observation suggested the possibility that noise might correlate better with a parameter 
in which the temperature rise was suppressed. Accordingly, a quantity called “APWEp^^j” 
was computed which simply omitted the 20 log AT term. Data for H6 were plotted against 
tliis parameter, as shown in Figure 21 . This figure makes clear that APWLpj^j is a poor 
predictor of measured noise compared to PWLpj^j. 

To determine if the scatter of the data of Figure 20 could be reduced, the poi its were 
segrepted into groups corresponding to runs at simulated idle, approach and takeoff con- 
ditions. Figure 22 shows that the idle data simply scatter about ±2 dB around a constant 
level as the burner exit temperature is increased by burning richer mixtures. Clearly PWEp^^^ 
fails completely to describe this behavior. 
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The approach condition data, Figure 23, have been divided into 3 groups as indicated by the 
broken lines. The three points on the line with slope of unity correspond to increasing fuel 
air mixtures, using only the secondary nozzle fuel system. Two points falling on an essentially 
horizontal line were obtained by supplying all fuel through the primary nozzle system. Com- 
pleting this picture, the vertical line connects two points that were obtained with both fuel 
systems operating; the lower point corresponding to supplying 35 percent of the fuel through 
the primary system, and the upper point representing a 50 percent fuel spbt. Taken together, 
these results indicate that the details of how fuel is supplied to the burner may have an im- 
pc'-tant effect on the noise. These details which presumably affect the time derivative of the 
heat release rate, should on physical grounds be of significance. By the same reasoning it 
cannot be expected that a simple prediction formula which is a function entirely of steady- 
state parameters can account for such variations. 

The observation that lower noise levels were observed at approach for configuration H6 by 
exclusive use of the second-^ry fuel system should not be generalized without having a deeper 
understanding of the unsteady burning process. For example the data in the cluster ai takeoff 
power, Figure 24, were recorded with a fuel distribution varying between about 10 and 35 
percent primary to total flow, and if plotted against percent primary flow indicate more 
scatter than a trend. Despite the need for caution, it b difficult to avoid spectulating that a 
special in-flight fuel system selection for use in terminal approach might be developed to offer 
a measure of combustion noise reduction capability. 

Tlie takeoff data, plotted in Figure 24, were taken in such a way that there was very little 
variation in temperature rise. Most of the points represent different fuel management pro- 
cedures. Therefore the takeoff data by themselves, as presented in Figure 24, do not pro- 
vide a good basis for evaluating 

Summarizing, for configuration H6, Eq. 9 of TMX 71627 is a reasonably good predictor of 
the overall variation with operating parameters of internally measured noise “power”. More 
detailed investigation showed that at some conditions the fuel management split produced 
noise variations that were not predicted by Eq. 9. In addition the acoustic “power” calcu- 
lated from the internal, near field pressure transducers is on the order of 16 dB liiglier than 
the predictions for far field power, a difference which is not unreasonable to expe< t between 
internal near field and far field quantities. 

Configuration S23 Results 

Results for tests on S23, which included runs at simulated idle, approach, climb, and takeoff 
conditions, may be described briefly as similar to those just described for configuration H6. 

A plot including all data of vs PWEp^^j, Figure 25, shows reasonably good correla- 

tion, the standard error of estimate being 2.4 dB as compared to a deviation of 3.1 dB from 
a horizontal line through the centroid. PWL^^^ averages about \ 2.6 higher than the far field 
PWL Suppression of the temperature rise e! ’ jct in predicting ^PWEp^^.^, , Figure 26, 
results in somewhat better correlation than obtained in H6, but offers no advantage over use 
of the more straightforward quantity PWLp^^j. 


Breakout of the data for idle, approach, etc. conditions is not as revealing as was the case for 
H6. In the idle data plot. Figure 27, the trend determined by the pair of right hand points is 
violated by the extreme left hand point. Examirution of the spectrum at this point, in com- 
parison with the idle cold flow spectrum indicates a very slight noise increase with burning, 
making calculation of relatively unreliable. This point, being dominated by aerody- 
namic noise should probably be discarded, in which case the two remaining points indicate 
good correlation between predicted and measured power. 

The approach condition data. Figure 28, do not span a sufficient range of either or 

PWL^ to establish information about conformity with or departures from predicted be- 
havior. Attempts to relate to fuel management, as was done for H6 are of question- 

able value due to the small range of PWL^^. 

Three points at different fuel- '.r ratios for the climb condition. Figure 29, show no variation 
as burner temperature rise is increased, although the change in the prediction parameter 
(about I.S dB) is probably too small to allow giving; much significance to this observation. 

At takeoff. Figure 30, the three data points spread out in both coordinates sufficiently well 
to allow observing that a slope 1 line would give a much better fit than a null effect or hori- 
zontal line. 

A general summary relating to test results for configuration S23 is that, as with H6, 
is a reasonably good predictor of the overall variation of noise with operating characteristics. 
The average difference between internally measured “power’* and predicted far field power 
is about 12.6 dB. Based on the data for the particular test points that were run, no informa- 
tion is available about the effects of other variables, sueb as fuel systems management. 

Results for Confin-.ations P-8, S-S, N-8, and N-JJ 

These configurations were tested in Phase 1 of the ECCP program at simulated idle and take- 
off condition only. Comments about the spectral characteristics of the internal pressure 
transducer signals are given in Reference 1 . 

Since the detailed testing and analysis performed for configurations H6 and S23 and reported 
immediately above have disclosed the general conformity of data trends with predictions 
and also possible departures therefrom, the current analysis of the Phase 1 data was limited 
to obtaining sample results. Specifically the four configurations listed above were selected 
from the Phase I program, and for each of these one idle point and one takeoff point were 
used in calculating power levels. The test points were selected to obtain comparable values 
of PWLpjjj among the configurations at idle and comparable values at takeoff conditions. 

Power levels for the above data are presented in Figure 31 . Despite the fact that only the two 
points for configuration P8 define a line having a slope close to unity, it should not be con- 
cluded that the other three configurations demonstrate systematic departures from applica- 
bility of PWL j as a prediction of variation of PWL^.^^ with operating parameters. It will 
be seen that the difference in PWL^^^j from idle to takeof*' ’ smaller (due to test procedures) 
than was the case for H6 and S23. Further, the plots for two configurations (Figures 
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20 and 25) show that due to data scatter, slopes very much dinerent from unity could be ob- 
tained by selecting just two points. 

The most significant information revealed in Rgure 31 is probably the general differences in 
mean levels among the configurations, an observation that was made in Ref. 1 on the basis of 
examining all the available spectra. These differences are discussed in the next subsection. 

Comparison Ammig Configurations 

From the relatively large number of data points available for cmifigurations H6 and S23 it has 
already been noted that the variations of with operating conditions is reasonably well 

predicted by has just been seen that insufficient data are presented for the other 

four configurations to either support or refute this observation. However, using the 2-point 
mean values for the Phase I configurations, as well as the more reliable mean values for H6 
and S23 it is possible to infer overall level differences among all the configurations. Table 11 
below gives the differences between corrected (observed) power level, PWL^ and predicted 
far field values, PWL^^: 


TABLE II 

MEASURED AND PREDICTED POWER LEVEL DIFFERENCES 


Configurations 

Difference: PWL^-PWL^^ 

H6 

16.2 

S23 

12.6 

P8 

23.3 

S8 

13.4 

N8 

10.2 

Nil 

6.3 

Average 

13.7 


Table II , -.ows that the level of PWL^^^^is significantly higher for all combustors than the far 
field noise predicted by referenced Equation (9). It is also seen that some combustors are 
appreciably quieter than others, assuming that PWL^^, derived from internal measurements, 
can be used to rank combustors for far field noise. Thus, while the variation of internal noise 
for a given burner was seen to be reasonably well described by the specific combination of 
operating parameters giving PWL^,^ , the levels of internal noise are appreciably higher than 
given by PWL^^ and there are large differences in level among various combustor designs. 
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These results are consistent with the previously discussed property of near fiekl measurements 
in producing high estimates of far fleld noise. Furthermore, while PWL ^ predicts the same 
noise for all combustors at a common operating condition, the measured results indicate large 
differences among configurations. These differences are hardly surprising, for a simple ex- 
pression involving only steady state parameters that does not account in some way for a variety 
of combustion chamber desigtu and fuel system details ignores the aero-thermodynamic de- 
tails that govern the unsteadiness of the heat release rate — the basic source of combustion 
noise. 

For background information a comparison of one-third octave band sound pressure ^ectra is 
presented in Figure 32. These spectra, which were used in conputing , are represen- 

tative of takeoff conditions for the six configurations of Table II and are the corrected spectra 
obtained by subtracting corresponding cold flow spectra from the takeoff point raw spectra. 

It may be noticed that the general decrease in spectrum levels conforms to the sequence of 
power levels given in Table II. An interesting feature is that all ^ectra contain signiflcant 
energy in the frequency range from about 400 Hz to 2000 Hz, which is appreciably higher 
than the range usually associated with combustion noise. A similar observation applies to 
the spectra at lower power conditions as may be inferred from the additional curves in Ap- 
pendix A-3 and those for other configurations tested in the Phase I program, contained in 
Reference 1. This frequency shift, which amounts to at least a full octave above previous 
data, cannot be attributed solely to general peculiarities of internally-measured pressure be- 
cause these previous data also include results of internal measurements. At this point no 
satisfactory explanation has been evolved to account for tliis result. 

A final pair of curves. Figures 33 and 34 completes the presentation of the pressure data. Al- 
though these figures pertain to data for configurations H-6 and S-23 their introduction has 
been deferred in order to facilitate presentation of results. These curves give for Configura- 
tions H-6 and S-23, respectively, plots versus PWL ^ of the uncorrected or raw power levels, 
denoted simply by PWL. It is evident from these figures that total “power”, PWL, uncor- 
rected for flow effects, does not correlate with PWLp^^j , since the data distribution is best 
described as scattering about an essentially constant value. Without giving statistical infor- 
mation about these distributions, the relatively poor correlation can be further appreciated 
by comparing these figures with the corresponding plots using corrected power, 

Figures 20 and 25. These last curves support the decision to avoid use of raw power level as 
a measure of combustion noise in this program. 

Internal Pressure Measurements and Far Field Combustion Noise 

Until a reasonable number of tests of a different type are conducted, the relation of internal 
pressure measurements to far field combustion noise will remain a matter of speculation. 

Even after correcting the data to reduce the influence of aerodynamic noise (using PWL 
a basic remaining problem is that transducers such as were used in this program, being 
positioned directly inside the source volume, measure the so-called hydrodynamic field rather 
than just a propagating acoustic fleld. In the area of fan and compressor noise, where the de- 
tails are comparatively well understood, it is known that pressure measurements taken very 
close to the rotating or stationary blades are meaningless as predictors of far field noise. On 
the other hand, measurements taken in the inlet and discharge ducts, at some distance from 
the blade rows can indeed be used to predict far field results (although the measurement pro- 
cedures and subsequent calculations are quite complicated). 
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These facts suggest that future combustion rig tests should employ transducer arrays in a dis 
charge duct sufficiently downstream of the zone of combustion so that the resulting measure- 
ments would give information about the propagating acoustic field. Such tests should also 
employ combustion chamber transducers, as were used in this program, to determine the de- 
gree of correlation between the combustion chamber measurements and the propagating com- 
bustion noise. (In the test arrangement used in the program reported here, it was not possible 
to install duct-mounted transducters in a location where reliable measurements could be ob- 
tained, due to the plumbing geometry and control valve locations employed in the combus- 
tion rig setup.) The objective of such tests would be to determine the degree of reliability of 
combustion chamber measurements systems in predicting far field noise, so that burner rig 
results could be used with confidence to develop combustion noise reductions, study the de- 
tailed mechanisms involved, and predict far field combustion noise for full-scale engines. 

5.2 TEMPERATURE INFORMATION 

As stated previously the objective of the fluctuating temperature investigation conducted in 
this program was to evaluate the potential for generating “indirect” combustion noise in a 
JT9D engine equipped with the test combustors, as a result of the interaction of the tempera- 
ture inhomogeneities in the combustor dischai^e stream with the pressure drop across a 
JT9D turbine. 

5.2.1 Method Of Analysis 

“Indirect” combustion noise, or “entropy” noise results from a mechanism whereby tempera- 
ture fluctuations or inhomogencities in a stream that encounters a strong pressure gradient 
produce acoustic waves. It may be helpful to regard tins mechanism in terms of local momen- 
tum fluctuations associated with the acceleration of fluid volumes of different densities im- 
bedded in the mean flow undergoing a change in velocity. Tliis mechanism was first identi- 
fied in the context of iircraft gas turbine noise by Candel (Ref. 5), although it is implied by 
the generalized unsteady flow analysis of Kovasznay and Chu (Ref. 6). Additional work dis- 
closing the importance of this source in core engine noise was performed and reported by 
Cumpsty and Marble, (ref. 7) and by Pickett (Ref. 8). The analysis used in this investigation 
assumes the pressure change to take the form of a sharp jump across an actuator disk, and 
was implemented by a computer program developed by Pickett based on Reference 8, which 
is enclosed as Appendix BI. 

The required inputs to the program, in addition to the steady-state turbine characteristics, are 
simply the mean square temperature fluctuations and the transverse and axial length scales 
of the tempciature distribution. These temperature field quantities were found for a total of 
1 7 operating points in the ECCP program by processing the signals from a fast response 
thermocouple system developed by Oils, References 9 and 10. 
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Tempentwre Field Characteristics 


Mean square temperature fluctuations were obtained from the thermocouple system voltages 
by means cf the conversion procedure described in the section on Acoustics Instrumentation. 
Although not required for input to the indirect noise calculation, 1/3 octave band thermo- 
couple voltage spectra were also obtained for general information and to obtain occasional 
chedrs of the overall mean square levels. In addition to the mean square intensity of tem- 
perature fluctuations, axial and transverse length scales are required in order to predict indi- 
rect combustion noise. The axial length scale, L^, was obtained from the autocorrelation 
function of the sampled thermocouple output voltages. In conformity with the mathemati- 
cal model of the indirect noise analysis, the characteristic ieng;. (or corresponding delay 
time) is defined as the point where the normalized correlation function has a value of I/e = 
0.37. Thus, for the axial length scale, the delay time corresponding to a normalized auto- 
correlation of 0. >7 was read from the autocorrelogram and multiplied by the stream ve.ccity 
to give L^. S^'vera! values of were determined from different thermocouples in the array 
at different rig operating conditions. 

Transverse length scales were obtained by cross-correlating paired outputs from transducers 
on 4 radial rakes. By cross correlating successive radial thermocouples with the end element 
of a rake and plotting the normalized cross correlation coefficients at zero time delay versus 
probe separation distance, the radial length scale was determined by entering the curve 
drawn through these points at a value of 1/e and reading out the corresponding distance. The 
transverse length scale could be determined by a similar procedure using corresponding ele- 
ments of different radial rakes. However, it was found early in the data reduction that the 
cross correlation between elements on the closest rakes was insignificant, due to the relatively 
large distance between rakes, as compared with the comparatively small separation between 
thermocouple elements on a common rake. Therefore the only inference that could be 
drawn was that the circumferential length scale was significantly less than the rake spacing 
(order of 6 to 7 cm). In the absence of further information, for purposes of calculation the 
circumferential and radial length scales were taken as equal. 

It should be remarked that the temperature data were first examined for evidence of a pos- 
sible fluctuating structure that would not produce indirect noise. It has been reported 
(Ref. 9) that some combustors operate with long streams of relatively cool air issuing from 
the chamber cooling wall holes, that wander randomly in the position where they penetrate 
the burner exit plane. Such motion will not produce indirect noise. Thermocouples in the 
exit plane would exhibit time-varying fluctuations as a result of such streams sweeping across 
the elements. If the inteipreter of the thermocouple outputs were unaware of the actual 
fluid motion involved he could impute the fluctuation to the passage of hot spots in a purely 
"xial direction, in the manner postulated for the indirect noise temperature input fleld, and 
.^uld proceed to use the inferred intensities and length scales to compute indirect noise as- 
sociated with a turbine pressure drop when, in fact, no such noise would be generated. 
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Accordingly, the tenq>erature data were examiiMd for the possible presence of temperature 
streaks sweeping transversely in random directions across the thermocouple array. Their 
presence would be indicated by cross correlation functions between pairs of thermocouples 
that peaked at delay times that would be related to the streak transvme velocity and the 
separation between thermocouples. If, for example such a streak were moving radially, the 
delay times for the peaks of the voltage cross correlations between an element at one end 
and successively distant elements would progressively increase or decrease, depending on the 
direction of motion. A similar succession of progressively delayed peak cross correlations 
between circumferentially di^laced thermocouples would iiulicate the presence of circum- 
ferentially moving streams of hot or cold fluid. 

On the other hand, if in fact the temperature non-uniformities were in the form of globules 
convecting along steady streamlines, the peak cross correlation between adjacent thermo- 
couples would occur at zero delay time, for the arrival of such globules at adjacent elements 
would be synchronous. In this test program, it turned out that such indeed was the case, so 
that the input temperature fields had the globular structure for which the indirect combus- 
tion noise model applies. 

Prediction of Indirect Combustion Noise 

The downstream indirect acoustic intensity flux was evaluated in each of the standard third 
octave bands for the passage of stream hot spots through the successive blade rows of the 
JT9D 2-stage high pressure turbine: first stator, first rotor, second stator and second rotor. 
Total intensity was obtained by adding the downstream contributions from each of the four 
rows, and no attenuation was applied for the multi-stage low pressure turbine or the trans- 
mission loss through the exhaust nozzle. The low pressure turbine generates insignificant in- 
direct noise due to the much lower pressure drops involved. 

Calculations were performed by IBM S/370-168, programmed to evaluate the triple integrals 
(over frequency and the two transverse wavenumbers) in Equation (8) of Appendix B-1. 

This equation, which is too lengthy to reproduce here gives the intensity between ^ecifled 
frequency limits which are taken as applicable to the standard 1/3 octave bands. The trans- 
verse wave numbers are involved because they enter into the spatial distribution of the acous- 
tic fields and account only for the portions of the energy that are in propagating modes, as 
distinct from decaying modes. 


Some appreciation of the significant factors governing indirect combustion noise may be dis- 
cerned if this equation is rewritten as follows: 



= const. (AP)2 d'^ f f 

j-'y2 

F dcj ' dk 2 dk y 

J, 7 , 
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In this abbreviated form it is readily seen that the acoustic intensity is proportional to the 
siuare of the pressure drop, AP, across a row, and also to the mean square value of the (non 
dimensional) entering stream temperature fluctuations, 0 root mean square (r.m.s.) fluc- 
tuating temperature divided by steady state average temperature. Thus, doubling either the 
stage pressure drop or the r.m.s. fluctuating temperature eadi produces a factor of four 
change in acoustic intensity or 6 dB in power level. 

The effects jf the other variables affecting noise are much harder to extract since these vari- 
ables are combined as complicated functions in both the integrand, F, of the above expres- 
sion and also in the limits of the two int^rations with respect to wave number. However, 
some insight is provided by considering the fecial case where the axial length scale is 
sufficiently long such that most of the resulting noise will be restricted to frequencies that 
are tow enough for only the plane wave mode to propagate. In this case, the intensity in a 
bandwidth 2Ao) ', centered at the frequency of peak intensity, u is given by Equation 12, 
Appendix B-1. Writing this equation in simplified form ^es: 

— ^ = const (AP)2 2 l* . 

2Ag,‘ Uto'„-A« ^ " 


This expression shows the previously demonstrated effects of pressure drop and temperature 
level (which are exact) and further reveals that the intensity is proportional to or 

to the transverse correlation area of the “hot spots”. 


The principal effect of axial length scale is to establish together with the axial stream velo- 
city, the frequency of peal; intensity, cjj^, which is on the order of the ratio of velocity to 
axial length scale. However L also enters into determination of intensity in several more 
complicated ways. Let equation (1 1) of Appendix B-1 for the intensity of a plane wave in 
a frequency band be simplified to the form: 



= const (AP)2 0^2 LyL'^ 



The overall intensity is obtained by taking cj 'j = 0 and Evaluation of the exponen- 

tial between these limits is just (minus) unity, giving 


I (overall) = const (AP)2 fl2 LyL'^ 
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It is now seen that the overall acoustic intensity in the plane wave mode varies inversely with 
axial length scale. The inverse nature of this dependence may be appreciated qualitatively 
by recognizing that as increases the momentum fluctuations across the blade row become 
progressively slower, thus diminishing the source strength. In the limit of the change 

in momentum across the row is constant (independent of time) and no sound is produced, 
which is evidently the case for steady hot or cold streaks of unlimited length. 

Going in the other direction, as decreases, the overall intensity increases by virtue of the 
above described 1 /L^ dependence. In addition, however, because the characteristic frequen- 
cies increase, acoustic energy begins to propagate in higher modes as well as in the plane 
wave. Mathematically, the more complex Equation 8 of Appendix B-1 must be employed, 
and it is not productive to try to trace the effect of decreasing beyond remarking that 

the double integral over wave number in Equation 8 reflects th'^ cumulative energy contri- 
butions of successively larger numbers of higher propagating modes. 

Thus, as axial length scale is reduced the overall effects are to increase both frequency of 
peak acoustic intensity and the overall intensity as well. There will then come a point where 
a basic assumption of the analysis begins to break down as length scale is reduced, because 
the analysis assumes existence of an ‘'actuator disk” where the pressure drop occurs over an 
inflnitisimal distance. This assumption is well suited to the case of shock waves or to a choked 
cascade, but is otherwise increasingly violated as the length scale of the hot spots becomes 
less than a substantial multiple of the axial extent of the blade row. If in fact the temperature 
axial length scale and the blade row thickness are comparable the calculation will over-pre- 
dict the indirect noise because the actual time for the hot spot momentum to change will be 
longer as a consequence of the finite distance over which the local momentum changes occur. 

The high pressure turbine stages (two stators and two rotors) of the JT9D operate unchoked 
and arc not much thinner than the axial length scales that were measured and will be described 
presently. Therefore it can be expected tliat the results of the analysis performed here will 
tend to ovcrprcdict indirect combustion noise in this installation. 

5.1.2 Results 

Applying the methods described above to the measured data the following results were obtained 
Measured Temperature Characteristics 

Fast response thermocouple information was obtained from ng tests on combustors H6 and 
S23 at a total of 1 7 operating points. These points are identified in the listing of rig perfor- 
mance data. Appendix A-1 . The majority of the temperature data were obtained for S23 be- 
cause this combustor is closely related to the “final” configuration evolved in the main F'x- 
perimcntal Clean Combustor Program and also because high thermocouple mortality severely 
limited the information available from the H6 combustor test program. 
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Although they are not used as speciHc input for ''alculation of indirect combustion noise, 

1/3 octave band thermocouple voltage spectra were obtained and are presented in Appendix 
B-2. Sample spectra for H6 and S23 are shown in Figures 35 and 36. It is worth noting that 
the spectral shapes are quite similar for all operating conditions. In the appendical plots oc- 
casional high frequency prominences are exhibited, but these are chargeable to spurious elec- 
trical system noise since the amplifier system response characteristic falls at a 12 dB per 
octave rate beyond 2000 Hz. 

The general peaking of the spectra in the frequency range from 500 to 1000 Hz and subsequent 
decline should be interpreted with some care because of the frequency response characteristic 
of the thermocouple and the compensating amplifier used in playback of the tape-recorded 
signals. As seen in Figure 1 5b the s>;tem response function Iregins to roll off at about 400 Hz 
and is 6 dB down at KXX) Hz. Therefore the true peaks of the temperature spectra will be at 
somewhat higher frequencies than indicated in the spectra presented here. Naturally, the 
system response characteristic also influences the length scale determinations which are ob- 
tained from correlation functions. The impact of this response characteristic on the length 
scales and on the calculated indirect combustion noise will be mentioned after the following 
description of the measured length scales and temperatures levels. It turns out that this high 
frequency roll off is of little practical significance because the calculated indirect noise spectra 
have peak frequencies that are already higher than normal for core engine noise. 

To obtain predictions of indirect combustion noise the temperature information requirements 
are the axial and transverse correlation length scales, and the mean square or root mean 
square (rms) levels of temperature fluctuation. 

The axial length scale is determined from the autocorrelation function of the thermocouple 
output voltage. Figures 37 and 38 present sample autocorrelograms at approach and takeoff 
conditions for configurations H6 and S23. These figures show the characteristic delay time 
corresponding to a value of the normalized autocorrelation of 0.37 (= 1 /e) and the subsequent 
translation to length scale by use of the stream velocity. A complete set of autocorrelograms 
is included in Appendix B-3. The effect of the system function high frequency roll off dis- 
cussed previously is to make the axial length scale obtained from the correlograms somewhat 
longer than the true value, since correlation delay time and frequency are inversely related. 

As previously described the transverse length scales, and , are obtained from plots against 
thermocouple separation distance of the normalized cioss correlation function at zero delay 
time for pairs of thermocouple output voltages. Figures 39 and 40 show for configurations 
H6 and S23 sample plots at approach and takeoff of the radial cross correlation, from which 
the radial length scales were determined. Since there are insufficient data points to justify 
fitting a Gaussian distribution function a simple decreasing curve drawn through the data is 
entered at an ordinate of I/e, determining the radial scale with sufficient accuracy. 



All the basic auto and cross-correlogratns are presented in Appendix B-3. The cross conelo- 
grams between adjacent pairs of circumferentially distributed thermocouples show normalized 
cross correlations on the order of less than O.OS, as may be seen from the appropriate curves 
in Appendix B-3. Since the closest circumferential spacing is on the order of seven times the 
radial thennocouple spacing (« 1 cm) and the radial cross correlation at 2 cm spacing is down 
to about O.OS, it is clear that the circumferential length scale is very much less than the closest 
circumferential spacing. In the absence of further information it was assumed as a matter of 
convenience that the circumferential and radial length scales were equal. Length scale infor- 
mation is summarized in the following tables, HI and IV. 

The flnal item of temperature information needed for indirect combustion noise is the mag- 
nitude of the temperature fluctuations. Root-mean-square values were obtained by true rms 
voltmeter during the correlation portion of the data reduction. A normalized temperature 
(rms fluctuation/steady mean absolute temperature) is used as a muLiplier to the output of 
the computer program, the input to which assumes for convenience a normalized value of 
one percent. Generally, several thermocouple voltages were read at each of the selected oper- 
ating points corresponding to simulated idle, approach, climb, and takeoff conditions. Tables 
III and IV summarize this information for configurations H6 and S23. Point weighted average 
values normalized to burner exit temperature are listed. Although not used for calculating 
noise the ratios of fluctuating temperature to temperature rise, AT, across the burner, are also 
tabulated for comparison with some previously acquired data. 


TABLE III 
COMBUSTOR H-6 

SUMMARY OF AVERAGE FLUCTUATING 
TEMPERATURE CHARACTERISTICS 



Length Scales 

Normalized RMS Temp. 

Condition 
Idle (2) 
Approach 
Climb (3) 
Takeoff 

Axial 

k 

cm (ft) 
5.2 (.171) 
4.9 (.159) 
4.7 (.155) 
4.7 (.154) 

Transverse 
L ,L 

y ’ z 
cm (ft) 

1.3 (.043) 

1.2 (.039 

1.0 (.034) 

1.0 (.033) 

% of Burner 
Exit Temp. 
d'c 

20.5% 

6.2 

3.9 

3.9 

% of Burner 
Temp. Rise 
(1) 

39.0% 

15.8 

6.9 

6.9 

(1) For reference only ; not used in noise prediction 

(2) Noise prediction not made for idle condition 

(3) Data not taken at this point. Interpolated values used for noise. 
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TABLE IV 
COMBUSTOR S-23 

SUMMARY OF AVERAGE FLUCTUATING TEMPERATURE CHARACTERISTICS 


Condition 

Axial 

K 

cm (ft) 

Idle (2) 

3.7 (.121) 

Approach 

4.45 (.146) 

Qimb (3) 

3.65 (.120) 

Takeoff 

3.9 (.128) 

Avg. of last 
3 pts (3) 

4.0 (.131) 


Transverse 


% of Burner 
Exit Temp, 
a'c 


11.7% 


% of Burner 
Temp. Rise 


21 . 0 % 


(1) For reference only; not used in noise prediction 

(2) Noise prediction not made for idle condition 

(3) Average of length scales at takeoff, climb, and approach used for noise predictions 


Commenting briefly on the temperature information in Tables 111 and IV, the following 
observations can be made: 


The axial and transverse temperature length scales are relatively insensitive to operating 
conditions and are generally similar in the two combustors measured. Measurement 
system response characteristics previously described have the effect of overestimating 
the true axial length scales. The effect of attenuation of high frequency information 
upon the transverse length scales is not so easily obtained. However, it is reasonable 
to expect that the true transverse length scales will either be unaffected or, as with the 
axial scale will be somewhat smaller than the indicated values. In any event, there is 
no basis to expect the transverse scales to be greater than the indicated values. 

On the other hand, the levels of the temperature fluctuations var>’ significantly with 
operating conditions, and this variation is quite different in the two combustors tested. 
For simplicity the idle data will be overlooked since this condition is of no practical 
concern from the noise standpoint. In configuration H6 the temperature fluctuations, 
normalized to burner exit temperature, drop from about 6% at approach to about 4% at 
takeoff. This change represents a reduction of about 4 dB in acoustic power. However, 
for S23 the normalized temperature rises from about 8% to 20% in this power change, 
which is equivalent to about a 7.5 dB noise increase. 
















3. The last column in the temperature tabulations expresses the rms fluctuation as a per- 
centage of the burner temperature rise in order to afford a comparison with some pre- 
viously acquired data. These data, presented in Reference 10, were acquired on com- 
bustors of three current production engines, and show fluctuations ranging from 7 to 
25 percent of the burner temperature rise. Apart from the idle condition, combustor 
H6 falls well inside the range of previously acquired data. On the other hand the 40% 
fluctuations measured from S23 at high powers are on the order of twice the values 
previously measured. These high levels have been ascribed to the relative proximity of 
the temperature measurement station in the burner exit plane to the combustion zone 
in S23, and have been judged to be consistent with more limited previous data taken up- 
stream of the discharge plane in conventional burners. 

4. In terms of indirect combustion noise power the temperature differences between H6 
and S23 amount to only about 2 dB at approach but at high power conditions S23 
should generate about 15 dB more noise than H6. The effects of the shorter axial length 
scales and larger transverse scales of S23 are such as to further increase its noise generat- 
ing characteristics relative to H6. 

Before discussing results of the noise calculation it should be mentioned that Appendix B-3, 
containing the temperature correlograms also includes several cross correlations (see pages 
181 and 182) between outputs of a sample thermocouple and a Kistler pressure transducer 
flush mounted upstream in the combustor S23. (See sections on pressure information.) 

These correlations were made to see if the internal pressur? fluctuations and discharge tem- 
perature fluctuations were causally related. The trivial values of the normalized cross corre- 
lation function of these quantities indicate that the time-variation of burner internal pressure 
and discharge temperature at the measurement locations are independent of each other. 

Predicted Indirect Combustion Noise 

Using the temperature information for combustors H6 and S23 summarized in Tables III and 
IV, the noise power generated by the interaction of the temperature inhomogeneities with the 
high pressure turbine of a JT9D engine was calculated from Equation 8 of Appendix B-1. As 
described previously, the downstream noise contributions of each of the 4 stages (2 stator 
and 2 rotor) were evaluated and summed. After passing through a stage the axial length scale 
was stretched in the ratio of exit to entrance axial velocities to establish the corresponding 
entrance length scale for the next stage in a manner that would conserve the . ‘me flux of the 
temperature eddies. 

Figure 41 presents the results of the indirect noise calculation for configuration H6 at approach, 
climb, and takeoff conditions listed in Tabic III. Downstream propagating power levels in 
1/3 octave bands are shown for each of the four high pressure tu'-binc blade (or vane) rows, 
together with the resulting spectral sum of these contributions. Similar results lor S23 are 
given in Figure 42. In both cases the length scales listed in Tables III and IV were used. How- 
ever, a common constant value of the normalized fluctuating temperature of one percent was 
used both as a matter of convenience and to illustrate the effects of variations in length scale 
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and turbine characteristics. Actual power levels corresponding to measured values of normal- 
ized fluctuating temperature, d'c, also listed in these tables, are obtained by merely adding 
20 log d'c to the one percent levels. On the basis of a common 1% temperature variation the 
overall power levels (integral of the resultant spectrum) of S23 range from about 3.S to * > 
dB higher than those of H6 as a result of its larger transverse and smaller axial length scales, 
in conformity with the background discussion in the pre /ious section. 

As seen in Figures 41 and 42 the second stator contributes on the order of 2 to 3 dB more 
than the next important stage, the first stator. The second and first rotors, in turn, are pro- 
gressively lesser contributors to the total power. These differences are a consequence of the 
’-tailed interstage variations of the turbine parameters entering into equation 8, Appendix 
B-1 and do not incorporate any assumptions for attenuation by downstream stages of noise 
generated upstream. In practice, of course, upstream noise will be progressively attenuated 
in its transmission through successive downstream stages, and the total resulting high pressure 
turbine noise also will be significantly attenuated by the 8 blade (or vane) stages in the low 
pressure turbine. Calculation of these attenuations and the further transmission loss at the 
exhaust nozzle is beyond the scope of the program, but rough estimates place the probable 
total attenuation between 5 and 10 dB. 

When the 20 log 0’c effect of the measured temperature fluctuations is added to the overall 
power levels obtained on a one percent basis the results listed in Table V iire obtained for 
combustors H6 and S23. 


TABLE V 


CALCULATED INDIRECT COMBUSTION NOISE 
JT9D ENGINE WITH CONFIGURATIONS H6 AND S23 
(OVERALL POWER LEVEL: dB re 10 WATTS) 



H6 

S23 

Approach 

129.6 

135.5 

Climb 

127.4 

145.8 

Takeoff 

! 27.7 

145.9 


The H6 noise levels decrease slightly with engine power as a result of the decline in levels of 
temperature fluctuations. For S23, on the other hand, the small increase with power on a 
constant one percent temperature basis is augmented by the large rise of input temperature 
fluctuation to net about 10 dB more noise at high power than at approach. 
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Evaluation of Indirect Noise Results 


Before evaluating the significance of the calculated levels of overall indirect noise some re- 
marks on the spectral distribution of this energy should be made. Examination of the spectra 
in Figures 41 and 42 shows that the indirect combustion noise power b a maximum in the 
range from 1000 to 1250 Hz, for both combustors over the power range from approach to 
takeoff. 

However, the class of phenomena called core engine noise, of which combustion noise is con- 
sidered an outstanding contributor, is characterized by spectra peaking well below 1000 Hz, 
usually in the 300 to 600 Hz range. This observed feature, common to several different types 
of aircraft turbine powerplants, is well documented in References 2 and 1 1 and in the liter- 
ature cited therein. 

For this reason, if the indirect combustion noise calc^dated for combustors H6 and S23 were 
to be actually generated in a JT9D engine undergoing acoustic tests, it would probably fail to 
be recognized as a core noise phenomena, assuming in the first place that its spectral levels 
were high enough to stand out from the broadband fan discharge noise spectrum. The cause 
of the higher than normal peak frequency of the predicted indirect noise is, of course, the 
relatively short axial temperature length scales measured for configurations H6 and S23. As 
in the case of the relatively high amplitudes of temperature fluctuation an explanation may 
be advanced that they are presumably a consequence of the detailed differences in combusto. 
design, with the relatively shorter than conventional Iciiiith of tiie com'-ustors probably being 
the most significant parameter. 

It may also be recalled that the burner internal fluctuating pressure measurements gave spectra 
(both as measured and when corrected for cold flow noise) having significant energy a :’' 
beyond the 300 to 600 Hz range. While these were characteristics of the vers' r. ' f . ;> .1 

it is not known how to transform them to obtain far field power spectra, it is of some interest 
to note that both the calculated indirect combustion noise and the measured internal pressure 
spectra differ significantly from t!'e conventional spectral representation of core engine noise. 

Finally, the probable significance of indirect combustion noise in the operation of a JT9D 
engine with the test combustors will be estimated by comparing the computed power levels 
with certain other core engine information items. 

Table VI gives farfield power levels of total low frequency core noise measured on several 
turbofan engines, including a JT9D, and also the values obtained by the recommended interim 
prediction formula. Equation 9 of Reference 2. Also listed are rounded-off values of the 
computed indirect combustion noise for burner configurations H6 and S23 operating in a 
JT9D engine. No low pressure stage turbine or exhaust nozzle attenuations have been applied 
to these two items, copied from Table V, but the measured and predicted levels do incorp- 
orate surh attenuation. 


TABLE VI 


APPROXIMATE TOTAL CORE NOISE 
POWER LEVELS 
PWLdBrelO-*^ WATTS 


item 

Approach 
Power Range 

High 

Power Range 

Engine Cl 

143 

155 

C2 

143 

(n.a.) 

C3 

140 

(n.a.) 

B 

140 

152 (est) 

A1 

144 

(n.a.) 

A2 

145 

(n.a.) 

Eq. 9, Ref. 2 

(turboshaft) 

(153) 

(162) 

turbofan 

145 

154 

Calculated, config. 



H6 

130 

128 

Calculated, config. 



S23 

136 

146 


The unavailable (n.a.) figures at high power operation for several of the above entries reflect 
the inability to detect or extract a spectral prominence in the expected frequency range. Levels 
from Equation 9 of Reference 2 for turboshaft engines are shown, together with correspond- 
ing values that are 8 dB lower, reportedly applicable to turbofan type engines. 

From Table VI it appears that a representative total core noise power level for turbofan en- 
gines at approach powers is about 143 oB and at high power is about 1 53. 
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Using these figures as representative standards for comparison it follows that: 

1 . The corresponding calculated levels for configuration H6 of 1 30 and 1 28 are at 
least 1 3 dB lower, so indirect combustion noise would not be significant for this 
burner. 

2. Calculated indirect combustion noise levels for configuration S23 are about 7 dB 
lower than representative total core noise. If only 3 dB low turbine and exhaust 
nozzle attentuation is assumed the indirect noise from S23 will be 10 dB less than 
total representative core noise. Therefore, for configuration S23 as well, indirect 
combustion noise is unlikely to be detected in the JT9D engine signature. 

Combining these evaluations of the overall levels of indirect combustion noise with the re- 
marks about the higher than normal peak frequency of their spectra it follows that indirect 
combustion noise is unlikely to be of practical significance for the JT9D and similar engines 
equipped with either of these experimental combustors. 
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& SUMMARY OF RESULTS 

!. Signals from transducers measuring pressure fluctuations inside combustion configura- 
tions were found to require special processing to provide a measure of combustion- 
generated noise. The overall sound pressure level indicated without such processing was 
significantly affected by sound generated aerodynamically. 

2. A measure of the near field of combustion noise was obtained by using only the portions 
of the transducer spectra in which combustion noise dominates. The lelevent spectral 
bands were found by comparing spectra for flows with and without combustion. In 
applications where cold flow cannot be produced in the combustor (as in an operating 
engine rather than a test rig), use of the spectra between about 100 Hz and an upper 
bound of 4K Hz would probably serve as a reasonable general guide. 

3. The resulting spectra were found to contain significant energy in frequency bands on the 
order of an octave higher than the spectra usually associated with combustion noise. 

This property applied essentially to all configurations and operating points. 

4. The sound pressure level determined as described above, and designated as OASPL^^, 
(overall sound pressure level, corrected) w^used to calculate a hypothetical combustion 
noise level (PWL^.^^^) using the relation P = p^ A/p c, where P, p, p, c, and A are values 

in the plane of the transducer of “Power” pressure, density, sound speed, and section 
area, when the local Mach number is relatively low. This “Power” is referred to as hy- 
pothetical because the equation from which it is computed implies that the acoustic 
field pressure and particle velocity are in phase, a situation that exists only in the “far 
field” (at large distances from the source). Due to neglect of this phase angle factor, 
calculations based on internal measurements overestimate the true acoustic power. 

5. For this reason, future experiments to measure combustion noise should provide pressure 
transducer stations downstream of the actual combustor, insofar as it is practicable to 

so locate them. In such cases, sufficient transducers should be employed to insure that 
higher duct modes as well as the plane wave are identifiable and taken into account 
when computing acoustic power. 

6. Extensive pressure data taken on configuration H6 and S23 were processed as described 
above to compute an acoustic power level called PWL^.^^, and these values were related 
to rig operating parameters by using Eq. 9 of TMX-71627. This equation predicts 
acoustic power level, called here PWEp^^j by the expression: 

PWL„,^ - 56.5 » 10 log,, m. (at ' 

where m^ is weight flow through the combustor, AT is combustion temperature rise, 

P&T are total Pressure and Temperature at the combustor entrance and the subscript 
amb refers to ambient or standard atmosphere values. 

A test of PWLpjjj , as a predictor of measured power was obtained by comparing 
“measured” acoustic power (PW^^^j, obtained from the transducers) with values deter- 
mined from Eq. 9 above. 
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7. It was found that the referenced Eq. 9 provides a reasonable prediction of acoustic power 
for Configurations H6 and S23, apart from the (xmstant. The measured power varies 
with mass flow, temperature rise, and combustor entrance pressure and temperature in 
the marmer described by referenced Eq. 9. However, measured power levels were ap- 
preciably higher than predicted, which is consistent with the above remarks that the 
near field transducer levels lead to overestimates of acoustic power radiated to the far 
Geld. 

8. Sample data at simulated idle and takeoff conditmns for 4 combustors tested in Phase 
I of this program, were insufficient to provide a valid judgement of the ability of Eq. 9 
to predict noise for these configurations. For the limited data used, measured power 
levels had a range of about 17 dB from the “loudes. o the most quiet, corrected to 
common rig operating parameters. 

9. Despite the general ability of to provide estimates of acoustic power variation, 

important deviations from predicted behavior were observed. For example, at idle con- 
ditions, tittle systematic change in measured noise was noted for some configurations as 
the temperature rise was increased by adding fuel. As another example, data taken on 
configuration H6 at approach conditions showed significantly different behavior as a 
function of temperature rise depending on how the fuel supply was proportioned be- 
tween primary and secondary nozzle systems. 

10. The deviations noted in the above paragraph are not anomalous since no simple predic- 
tion formula that fails to provide inputs in the form of aero-thermodynamic field micro- 
structure can account adequately for combustion noise. One of the outstanding source 
terms of such noise is the time rate of change of heat release rate. This quantity is cer- 
tainly affected by details of how the fuel is injected and burned, and the turbulence 
structure of the flow, and no provisions exist in any simple formula for quantifying 
these crucial parameters. 

1 1 . Using fast response thermocouples in the burner discharge stream, fluctuating tempera- 
ture characteristics were determined for combustors H6 and S23 in order to compute 
indirect combustion noise. The required temperature parameters were axial, radial, and 
circumferential length scales of the temperature inhomogeneities and the root mean 
square (rms) amplitude of the fluctuation. In addition, sample one-third octave tem- 
perature spectra and a limited number of cross correlations between a thermocouple 
and an internal pressure transducer were obtained. 

1 2. For both configurations the temperature length scales remained essentially constant over 
the operating range from simulated idle to takeoff. Also, corresponding length scales of 
both combustor streams were comparable. The axial length scale was about S cm for 
H6 and about 4 cm for S23. Radial transverse length scales were approximately 1 cm 
and l.S cm for H6 and S23 respectively. Circumferential length scale could net be 
directly measured since the spacing between adjacent thermocouples was substantially 
larger in the circumferential direction than radially; it was taken equal to the radial 
length scale for input to the noise calculations. 
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1 3. Thermocouple voltage spectra are generally similar in shape for both combustors over 
the operating range. However the maximum one third octave band levels occur at freq- 
uencies that are on the order of an octave higher than are usually associated with com- 
bustion noise or low-frequency core engine noise. This characteristic is reflected in the 
magnitude of the axial length scales and subsequently in the spectra of the computed 
indirect combustion noise. 

14. To determine if the discharge temperature fluctuations were causally related to the 
combustion chamber pressure variations, sample cross correlations were made between 
signals from a hot wire and from one of the flush-mounted pressure transducers. Es- 
sentially zero normalized correlation coefficients were noted at all sampled operating 
points, indicating that the temperature and pressure fluctuations were completely inde- 
pendent processes. 

15. To facilitate discussion of the root mean square (rms) amplitude of the temperature 
fluctuations, this amplitude is expressed as a percentage of the steady-state burner tem- 
perature rise. Excluding the idle operating point, which has little significance from the 
noise standpoint, combustor H6 temperature levels decreased from about 16% at ap- 
proach power to about 7% of the burner temperature rise at high power operation. 

These figures are well within the 5% to 25% range previously measured for three different 
production engine burners. Configuration S23 behaved quite differently, increasing 
from about 1 7% temperature fluctuation at approach to 41% at takeoff power. Tims 

at low power (approach) H6 and S23 develop comparable temperature fluctuations that 
fall within the range of other, production type burners. At high powers, however, the 
temperature fluctuations of S23 are significantly higher than those of H6 and other 
production engine burners. 

16. Calculations were made to estimate the indirect combustion noise that would result 
from the passage of the temperature fluctuations described above tlirough the higli pres- 
sure turbine of a JT9D engine. Tlie shapes of the resulting acoustic power spectra for 
combustors H6 and S23 were quite simitar and did not change significantly from ap- 
proach to takeoff power conditions. The one third octave band spectra peaked in the 
range of 1000 Hz to 1 250 Hz, which is on the order of an octave higher than spectra 
usually associated with combustion noise or low frequency core engine noise. 

17. Levels of calculated indirect acoustic power, however, differed significantly between 
the two combustors and also varied with operating conditions in different ways. For 
H6 the overall acoustic power level was about 130 dB (re 10'*^ watts) at approach and 
dropped slightly to 1 28 dB at takeoff power, primarily as a consequence of the reduced 
temperature fluctuation. In configuration S23 the calculated indirect noise power 
level at approach was about 136 dB and increased to 146 dB at high power, reflecting 
the increase in temperature fluctuation. 
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1 8. Significance of these levels was evaluated by comparison with previously determined i 

total engine low frequency core noise obtained from measurements on turbofan engines | 

and from the interim prediction procedure discussed in the section of this report on ' i 

internal pressure measurements. Representative total core noise power levels were taken | 

as 143 dB for approach and 1 S3 dB at higli power. I 

If the noisier of the two configurations, S23, is considered, and if it is assumed that only 

3 dB of attenuation of indirect noise is provided by transmission through the multi-stage . > 

low pressure turbine and the exhaust nozzle, farfield indirect noise power levels of 1 33 

dB at approach and 143 dB are projected. ^ 

< 

19. Comparing these projected levels for S23 with representative figures for turbofan engine 

low-frequency core noise it is seen that they are about 10 dB lower than total core noise \ 

levels from all sources of typical turbofan engines. Levels for the other combustor, H6, 

will be significantly lower yet. Since the input measured temperature fluctuations for 
S23 were appreciably higher than indicated by previous experience, it appears that in- 
direct combustion noise for any of the experimental clean combustors can be ruled out 
as a significant component of JT9D engine core noise. 


7. CONCLUSIONS 

1. Pressure fluctuations inside a combustor, when corrected for the effects of aerodynamic 
noise, may be used to compute a rough estimate of acoustic power due to unsteady 
burning. The effects of near field transducer location and neglect of transverse mode 
structure are such as to lead to an overestimate of the true acoustic power. 

2. For a given combustor, the changes in internally measured acoustic power with varia- 
tions in burner operating parameters are reasonably well predicted by Equation 9 of 
NASA TMX 71627. 

3. For the combustors tested, the average of the internally measured power is about 14 dB 
more than predicted by this equation for the far field, a difference which is consistent 
with the complications of internal near field measurements noted in item 1 above and 
neglect of turbine and nozzle transmission losses. 

4. Differences in noise level of up to 15 dB were observed among different combustors, 
indicating that design details not specified in a simple expression such as referenced 
Equation 9 have important effects. Further, in a particular burner design, the effects of 
such variables as fuel supply system management may cause the variation of noise with 
operating parameters to depart significantly from the predicted variation. 
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5. Measurements indicate that the Swirl Vorbix Combustor, S23, produces discharge tem- 
perature fluctuations that are significantly higher than those previously observed on 
several production burners. Tlie Hybrid Combustor H6 produced lower levels of fluctua- 
tion that were more in line with previous experience. 

6. For both combustors, however, indirect combustion noise produced by interaction of 
the temperature fluctuations with the high pressure stages of a JT9D turbine was cal- 
culated to be significantly (10 dB) lower than observed values of total low-frequency 
core engine noise from all sources. 

7. Since these low sound power levels included the projection for an experimental burner 
with higher than normal discharge temperature fluctuations, and also since the resulting 
noise spectra peaked more than an octave above typical core noise spectra, it is unlikely 
that indirect combustion noise could be detected in the operation of engines such as 
the JT9D. 
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Figure la Swirl Can Combustor - Baseline Configuration (N-SeriesJ 
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Figure Ic Swirl Can Combustor - Combustor Modifications (Ref. Figure la) 
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Figure 4a Swirl Forbix Combustor S-11 (Phase II) 
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Figure 6b Uxterior View of Basic Hybrid Combustor 
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Figure 7 Hybrid Combustor Configuration H-6 (Phase H) 
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Figure 8a Schematic of Combustor Rig 
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Figure 10 Typical Transducer Installation (detail) 




















VOLTAGE RATIO 



Figure 15a Schematic Diagram of Thermocouple System Circuit 



Figure 15b Overall Frequency Response of Thermocouple and Compensating Amplifier 
System (Also Showing Sample Uncompensated Thermocouple Response) 
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Temperature Axial Length Scales From Auto Correlograms : S-23 
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APPENDIX A-1 


RIG PERFORMANCE AND POWER LEVEL DATA 

This Appendix contains tabulations of rig performance para- 
meters from Phases I and II of the Experimental Gean Com- 
bustor Program, and identifies configurations and opr rating 
point numbers. Acoustic power levels are also tabulated. 



RIG PERFORMANCE DATA > REF 1 
List of Symbolt 


SUBSCRIPTS 


4 

5 
t 
s 
i 
b 

pri 

sec 

SYMBOLS 


P 

T 

Wa 

Wf 



station upstream of burner. 

station downstream of burner 

total condition 

static condition 

ideal 

burner 

primaiy zone 

secondary zone 


Pressure 
jmperature 
mass flow rate of air 
mass flow rate of fuel 
fuel-air ratio 

mass flow parameter 


UNITS 

ATM 

“K 

kg/sec 

kg/sec 

kgy/i^ 
sec - m* - ATM 


NOTES: 


- 

M = 

A4?t4 


, where A 4 is defined as 0.02002 m^ for a 90 


sector JT9D burner. 
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k computed from carbon balance Aiel-air ratio, corrected for unrcacted species (CO. THO present in exhaust {*>• 



TEST RESULTS FOR HYBRID COMBUSTOR 
CONFIGURATION H3 





















































TEST RESULTS FOR VORBIX COMBUSTOR 
CONFIGURATION f23 
FROM REF 4 
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TABULATIONS FOR POWER LEVEL - 360° ANNULUS 
(UNITS AS SPECIFIED BY NASA TMX - 71627) 



PO'>it 


We (90°) 

"6^ 



”*‘■”'•'‘.•13 
dBre 10^-* 

cor 

PWL 

Cnnfig. 

No. 

Condition 

Ibm/tec 

°R 

V’o 

V^O 

watts 


N-8 

1 

Cold Idlf 

13.07 

0 

3.10 

1.50 



143.8 


32 

Idle 

13.03 

793 

3.77 

1.58 

13" .2 

153.6 

154 .8 


0 

Cold T/0 

15.21 

0 

6J0 

2.67 



154.4 


79 

T/0 

15.28 

1167 

6.93 

2.66 

144.0 

150.8 

156.2 

N-11 

1 

Cold ldl« 

13.40 

0 

3.74 

1.58 



151.9 


33 

Idle 

13.38 

697 

3.75 

1.57 

138.2 

148.3 

154.4 


60 

Cold T/0 

14.02 

0 

6.19 

2.67 



149.5 


58 

T/O 

13.93 

1393 

6.16 

2.67 

144.1 

146.9 

153.3 

P8 

0 

Cold Idle 

9.61 

0 

3.74 

■>.66 



153.1 


11 

Idle 

9.61 

824 

3£9 

1.59 

137.9 

160.7 

162.5 


102 

Cold T/O 

15.30 

0 

6.78 

2.67 



158.9 


410 

T/O 

14.42 

10S6 

6.81 

2.68 

143.0 

166.8 

167.6 

S8 

1 

Cold Idle 

12.01 

0 

3.53 

1.08 



149.3 


21 

Idle 

11.66 

823 

3.73 

1.56 

139.0 

154.5 

156.0 


2222 

Cold T/O 

15.43 

0 

6.93 

2.68 



149.9 


810 

T/O 

14.84 

1319 

6.84 

2.60 

145.0 

156.5 

■:58.9 

H6 

1 

Cold Idle 

8.75 

0 

2.92/ 

1.484 



158.0 


6 

Idle 

8.67 

1063 

2.927 

1.184 

138.3 

153.5 

161.2 


5 

idle 

8.72 

854 

2.927 

1.484 

136.4 

155.6 

160.2 


7 

Idle 

8.83 

550 

2.927 

1.484 

132.7 

154.2 

159.9 


8 

Idle 

8.43 

723 

2.927 

1.484 

134.8 

151.5 

160.1 


2 

Cold App 

17.50 

0 

6.807 

2.060 



155.4 


12 

App 

17.61 

681 

6.807 

2.050 

142,1 

159.3 

161.9 


11 

App 

17.74 

490 

6.807 

2 .osr 

' 0,2 

158.9 

161 3 


15 

App 

17.43 

700 

6 807 

2.053 

■..!2.2 

155.5 

159.7 


14 

App 

17.52 

490 

6.807 

2.050 

139.2 

153.5 

158.5 

13 or 

131 

App 

17.41 

365 

6.807 

2.050 

136.6 

150.4 

157.6 


17 

App 

17.41 

920 

6.807 

2.050 

144.6 

159.7 

163.3 


18 

App 

17.68 

881 

6.807 

2.060 

144.3 

155.0 

160.2 


4 

Cold T/0 

15.12 

0 

6.307 

2.667 



153.9 

23 or 

221 

T/O 

15.02 

1370 

6 307 

2.667 

145.1 

160.5 

164.0 


24 

T/O 

15.13 

1006 

6.807 

2.667 

142.5 

15“./ 

IBu.O 


24 

T/O 

15,13 

1006 

6.807 

2.667 

142. p 

159.7 

162 9 


19 

T/O 

15.28 

1347 

6.807 

2.667 

145.1 

161.6 

164.t 


20 

T/O 

15.11 

1373 

6.807 

2.667 

145.2 

161 6 

163.9 


21 

T/O 

15.15 

1355 

6.807 

2.667 

145.1 

loi.C 

16:l.O 


22 

T/O 

15.20 

1359 

6.807 

2.61. 

145.1 

162.7 

1-1.3 
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Config. 

S23 


TABULATIONS FOR POWER LEVEL - 360° ANNULUS 
(UNITS AS SJ^ECIFIED BY NASA TMX - 71627) 


Point 

No. 

Condition 

Wa(90 ) 
Ibm/ftc 

°R 

P^P 

4 o 


'^'-pr- 
dBr* 10’’^ 

cor 

Watts 

PWL 

1 

Cold Idio 

8.60 

0 

2.927 

1.484 



153.6 

2 

Idk 

8.92 

1059 

2.927 

1.484 

138.4 

152.6 

158.4 

3 

• 1* 

a70 

868 

2.927 

1.484 

136.6 

150.2 

156.1 

4 

It . 

a57 

405 

2.927 

1.484 

130.0 

148.9 

154.8 

181 

Cold App 

17.58 

0 

6.807 

2.050 



152.1 

5 

App 

17.66 

1013 

6.807 

2.050 

145 .p 

157.0 

159.7 

6 

App 

17.62 

918 

6.807 

2.060 

144.6 

156.9 

159.4 

7 

App 

17.68 

795 

6.807 

2.060 

143.4 

155.3 

158.8 

8 

App 

17.60 

1014 

6.807 

2.060 

145.5 

157.5 

159.5 

9 

App 

17.72 

953 

6.807 

2.050 

145.0 

156.7 

160.3 

10 

App 

17.77 

961 

6.807 

2.050 

145.1 

156.5 

159.2 

17 

Colo Climb 

15.68 

0 

6.807 

2.563 



152.7 

14 

Climb 

15.83 

1325 

6.807 

2.563 

145.4 

157.5 

158.3 

15 

Climb 

15.81 

1265 

6.807 

2.563 

145.0 

157.5 

158.3 

16 

Climb 

15.79 

998 

6.807 

2.563 

143.0 

157.2 

158.3 

175 

Cold T/O 

15.29 

0 

6.807 

2.667 



153.3 

13 

T/O 

15.17 

1358 

6.807 

2.667 

145.1 

156.2 

158.9 

12 

T/O 

15.12 

997 

6.807 

2.667 

142.4 

153.2 

158.6 

11 

T/O 

15.08 

900 

6.807 

2.667 

141.5 

151.1 

156.7 
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APPENDIX A-2 


TABULATED PRESSURE SPECTRA 

Pressure transducer one-third octave band spectrum levels 
are ta'^ulated in this Appendix for configurations H-6 and S-23. 



TABULATED 1/3 OCTAVE BAND 
TRANSDUCER PRESSURE LEVELS 
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TABULATED 1/3 OCTAVE BAND 
TRANSDUCER PRESSURE LEVELS 
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TABULATED 1/3 OCTAVE BAND 
TRANSDUCER PRESSURE LEVELS 
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ftEPfU)Dt]CffitLnY OF THE 
GltlOmA^ AGE IS POOR 


> 


TABULATED 1/3 OCTAVE BAND 
TRANSDUCER PRESSURE LEVELS 
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APPENDIX A-3 

PRESSURE SPECTRA PLOTS 

Machine plots of one-third octave band q;>ectnmi levels for 
H-6 and S-23 pressure transdw:er signals are contained in this 
^pendix, together with selected narrowband spectra. 
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APPENDIX \A 


SAMPLE CORRECTED SPECTRA 

Curves are presented here showing recorded pressure spectra, 
corresponding cold flow spectra, and the resulting corrected 
spectra used in calculating corrected acoustic powerlevels, 
PWLj,q|.. The method of correction is described in the main 
text section. Results and Discussion — Pressure Information - 
Method of Analysis. 
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APPENDIX B-1 


INDIRECT COMBUSTION NOISE THEORY 

The following material which concerns the analysis of indirect 
c(Hnbustion noise is the basis of the computer program used in 
the ECX:P noise program. This material is an edited version of 
the paper, “Turbine Noise Due to Turbulence and Temperature 
Fluctuations”, by G. F. Pickett, P&WA Acoustics Group 
(Reference 7). The paper was presented at the Ei^th Inter- 
national Congress on Acoustics, London, England, 23-31 
July 1974. 



TURBINE NOISE DUE TO TURBULENCE 
AND TEMPERATURE FLUCTUATIONS 


Abstract 

In addition to turbomachinery noise generated due to high 
speed unsteady flows interacting witli blades and vanes, signi- 
ficant low frequency noise can be generated by the passage of 
velocity and temperature “eddies” through large pressure gra- 
dients such as exist across turbine blade and vane rows. This 
noise generating mechanism may contribute to core engine 
noise observed from high bypass ratio turbofan engines. 

An analysis is presented that determines upstream and down- 
stream radiated noise generated by vonicity (turbulence) and 
entropy (temperature) waves interacting with a blade row si- 
mulated by an actuator disc. Radiated noise intensity is found 
to depend on the change in pressure across the rotor, and the 
intensity and transverse correlation lengtli scales of the tem- 
perature and axial velocity fluctuations. Axial correlation 
length scales together with the axial convection velocity deter- 
mine the frequency of noise generated. 


Introduction 

As fan noise from high bypass ratio turbofan engines is progressively reduced with in- 
creased inlet and aft duct treatment, and as jet noise is minimized by appropriate choice of 
engine cycle and bypass ratio, low frequency core engine noise and high frequency turbine 
noise become importan! sources of engine noise. Tliere is now substantial evidence (1,2) 
showing that noise from low velocity jets cannot be accounted for by jet noise .ilone, and 
that there exist sources of noise upstream of the jet nozzle which contribute significant lev- 
els of noise at frequencies generally similar to those produced by jet noise. 

Various mechanisms for low frequency core engine noise have been proposed, but, as 
yet, no definitive tests have been conducted that directly relate core engine noise to one spe- 
cific dominating mechanism. Following the example of fan noise, it is reasonable to suspect 
that unsteady flow velocities interacting with blades, vanes, struts and nozzle exit are signi- 
ficant contributors to the overall noise. An additional mechanism exists, however, that 
involves the large fluctuations in temperature and axial velocity generated in the burner pass- 
ing through the very large changes in pressure and swirl velocity across successive turbine 
stages. Physically, acoustic pressure waves and vorticity waves can be generated by this 
mechanism because axial velocity fluctuations and variations of density due to temperature 
fluctuations convecting with the fluid, experience abrupt changes as these fluctuations pass 
through large pressure gradients such as exist across turbine blade rows. 



In this paper, an analysis is presented that involves the solution of the three-dimensional 
linearized equations of continuity, momentum, energy, state and second law of thermodyna- 
mics. These equations are solved in conjunction with the satisfaction of appropriate steady 
and unsteady jump conditions across a vane or blade row that is modeled as an actuator disc. 

The general analysis includes noise generation due to fluctuating lift and drag forces at 
the disc in addition to the noise generated by fluctuating temperatures and velocities convect- 
ing through the disc. 

Total noise estimates from all the blade rows c. the turbine that contribute significantly 
are calculated by assuming that no coupling exists between the bl.ide and v;>ne ro.vs Under 
certain conditions Cumpsty and Marble (3) found that the coupling could be important and 
the present model will need to be extended to include this effect. An additi-aaal effect that 
could sigf’ificantly change the frequency distribution and sound power levels from any core 
engine noise mechanism is the duct exit impedance. Tliis impedance has been shown '4.5) to 
be large for low frequencies and to change with flow speed through the duct exit. The ef- 
fect of the detailed flow at the tailpipe lip on the impedance has not. to the author’s know- 
ledge, been considered in any analysis, but it probably does change the impedance at the 
duct exit. If this 'itter effect is significant, it should be expected that factors such as for- 
ward speed, fan duct geometrv and split or common flow nozzle, will affect levels of core 
engine noise. 

ANALYSIS 

Analytical Model 

The method considers small unsteady perturbations of the mean flow through an ‘‘ac- 
tuator di.se” blade or vane row that represents cither a turbine rotor or nozzle guide vane 
stage. The analysis also accounts for the large changes in mean scream flow variables that 
are not only an inherent design feature of turbine stages, but also a key factor in the noise 
generating mechanism considered in this paper. 

Consider a cascade as shown in Figure 1 , where the turbulent velocity superimposed on 
the mainstream flow (x, t) entering the cascade is u^^ (x, t) and where temperature fluc- 
tuations (x. t) convecting with the mean stream are superimposed on the upstream am- 
bient temperature Oo, • In the Appendix it is shown how these perturbations convecting 
through an actuator disc representation of the cascade, together with fluctuating lift and 
drag forces at the cascade give nsc to acoustic waves. Neglecting the noise due to unsteady 
lift and drag at the cascade, which has been considered in a .separate paper (6), the expression 
for upstream and downstream acoustic pressure waves due only to temperature and axial 
velocity fluctuations traversing the cascade reduces to; 
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Pp^Cu.k^.k,) = (w.k^.k^) + 


u. 


p„ U k <'> 

*^0 X 


»2 


APo ■*■ ^ A V 


Uc^(w.k^.k,) 


(I) 


A A 


where i = 1 and 2 represent upstream and downstream respectively and where the acoustic 
waves and the temperature and inflow turbulence are related in frequency and wave num- 
ber space. The values (i, j = 1 2) and AA are written out in full in the Supplement of 
this Appendix and other terms are as described in the Notation. 


INFLOW VELOCITY 




Z ~ SPANWISE 


figure 1 Cavciii/c Representation of Turbine lilade or lane Row 


PAUE IS POOR 
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The upstream and downstream acoustic pressure spectra can be calculated from equa- 
tion (1) by forming the product of the acoustic pressure wave in frequency /wave number 
space with its appropriate complex conjugate p* (co, , k^). The resulting- spectra will be 
seen to depend not only on the spectra of the inflowing temperature (o>, k^ , k^ ) P and 
axial velocity fluctuations ju^ (cj, k^ , k^ ) 1^ , but also on terms coupling tliese fluctuations, 

6 ^ u*^^ and 0 * u^ . Once the pressure spectra are known, tlie upstream and downstream 
acoustic intensity Radiating per unit area can be calculated from equations (A30) and (A31 ) 
in the Appendix by summing over tnose wave numbers which give rise to propagating ac- 
oustic modes. 

From the form of equation (1 ), it can be seen that noise is generated by temperature 
fluctuations only when there is a non zero pressure change across the actuator disc and by 
axial velocity fluctuations when there are non zero changes in either density, mean swirl 
velocity or mean axial velocity across the disc. In order to derive the complete noise field 
due to these mechanisms, the fluctuating axial velocity and temperature spectra in frequency/ 
wave number space and the cross spectra between axial velocity and temperature fluctuations 
are required. Because of the difficulty involved in providing appropriate inputs, only the 
noise due to convecting temperature fluctuations traversing the actuator disc will be consi- 
dered in this paper. Although the noise due to other terms may not be negligible, they 
would exhibit results having properties similar to those that will be brought out by consider- 
ing only the noise due to convecting temperature fluctuations. In most cases, this mechanism 
probably dominates that due to axial velocity fluctuations because, in general, turbine stages 
are designed with very large pressure drops and changes in mean swirl velocity, but only rel- 
atively small density and mean axial velocity changes which make the first two terms in the 
square bracket of equation (1) small. Tlie term containing the large swirl velocity term, how- 
ever, is multiplied by the circumferential wave number k^, . Since the mechanism be; ig con- 
sidered is essentially a low frequency noise generator, only the plane wave mode (i.e. k^ = 0) 
propagates over the lower end of the frequency range of interest in low frequence core noise. 
In this range, no contribution can be made by the term involving large swirl velocity changes. 

Considering, therefore, only noise due to convecting temperature fluctuations, the up- 
stream and downstream pressure spectra become. 



Thus, for a specified frequency/wave number spectrum of the convecting temperature fluc- 
tuations at the actuator disc, the upstream and downstream acoustic intensity flux due to 
these fluctuations traversing the disc can be calculated. 
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Representation of Fluctuating Temperature Input 

In order to make the infinite two-dimensional actuator disc of Figure (I ) be a more re- 
alistic representation of a turbine stage, the disc is taken o rcfiect the azimuthal 2 rr radian 
periodicity inlierent in actual rotor or stator stages. The temperature distortion entering the 
actuator disc representation of a turbine blade row is thus taken to be repetitive over a period 
of one circumference. 

In addition, the temperature distortion is assumed “frozen” as it convects through the 
actuator disc and to have a cross correlation in a coordinate system moving with the free 
stream of the form 


(x y z ) = 01 exp - 



( 3 ) 


where (x', y', z') represents the separation distance between two reference points. and 

Lj are representative of the axial, circumfc^ntial and span wise correlation length scales of 
the convected temperature distortion and 0 ^ is the mean square intensity. 

Applying the condition of periodicity over the range 2irr (where r is representative of 
the radius of the turbine stage) and assuming that L^, is much smaller than one circumfer- 
ence, the cross correlation of the temperature distortion can be written as a Fourier series 
in y' as 

ij y’ 

oo 

Rq (y', z', t ) = ,1 A. (z', t) e ^ , 

where 

■ . I 

y _ tjy 
L' f 

e e dy ; 


and x' has been replaced by U, r. 
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Since < < 2jrr. the limits of integration .an be extended to infinity permitting the above 
to be evaluated as 


U? 


i’ L’ 

J S ij 


ir'' L « 1 2 12 

R(y'.4r) = -— ^ 0^ e c 

2irr ' 


£ e 


4r* f 

e (4) 


The spectrum of the inflow temperature distortion is then readily calculated by Fourier 
transforming in space and time, so that 


I® (w, ky.k,)|* = 


■ L* 1 


8i’/*U,r 


1-- , if , 


(-. - i)- 


From the form of this equation, the following non-dimensional parameters (denoted 
by primes) are suggested; 


v-^ 

S r 


h T 


y = . 


If = r k 


k; = rk, e: = 


Thus, in non-dimensional parameters, 


|g(co.k,.k,)i^ = < H«.,'.k;,k;) 


where 


H(w'.k;,k;) = l; l; l; e 




L', 


«(k;-j) 
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Results and Discussion 


The overall sound power radiated to the far field due to convecting temperature fluctua- 
tions interacting with turbine blade rows can be obtained by substituting equations (6, 7) for 
the temperature distortion in equation (2) to obtain the acoustic pressure spectrum and then 
applying equations (A30) and (A31)of the Appendix to determine the acoustic intensity 

flux. 


If the radiated intensity integrated over the non-dimensional frequency band [cj’ , ca' 1 
b written in terms of non-dimensional parameters, equation (A3 1) for the upstream (i = 1 ) 
and downstream (i = 2) acoustic intensity flux becomes: 


kO)' k<*>' 

Iti tmi 



(AP)> p 

* 

* ii 

/ /. 

t HR, B*.dw'dk;dk;, (8) 

(O 

1 


to; 

•{(1)' 



1 *1 


where 




# 

r -Wy + 

(■ly 1 


kO) = 

M tj' li— 

*1 



‘ »- 

M* -M^ 



L 

*1 •'I J 


k«)' = 

w' 

*i 



*1 

>/l -M' -M' 

*i >^1 






In this equation, the triple integral is non-dimensional since the terms retained in the integral 
are the non-dimensional input temperature distortion H (w k ^ , k \). the non-dimensional 
radiation factor R. (w k^ , k '^ ) and the term B. ^ (oa k^ , k'^ ), which relates the temperature 
H>ectrum to the acoustic pressure spectrum. The non-dimensional limits of integration over 
the wave numbers are obtained from the Appendix and are such as to only account for the 
propagating acoustic modes that are generated. 

In order to determine the general properties of this result the triple integral is evaluated 
numerically to account for a wide variety of operating conditions. It is constructive, however, 
to consider first a particular approxim.ate case, valid at low frequencies, that can be written 
out explicitly. 
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Since the temperatures in the turbine are very high, the speed of sound is as much as 
twice that of the ambient speed of sound and tire flow is characterized by relatively low 
Mach numbers. Tire wave number limits of integration in such cases reduce approximately 
to ± w* M« j. Now, if only very low frequencies are considered, so that I Max (k^ i = 

I it)* M»,| < 1, only the circumferential mode k'^ = 0 can propagate and only the j = 0 term 
in equation (7) contributes to the total. Also, at low frequencies there is only a small range 
of k *s that give rise to propagating modes (i.e., | I < 1 tj' M,, I). Over this range, H. It 
and vary only slightly and equation (8) simplifies including the expression for H from 
(7) to: 



<r.xb? 3 >l'.l;l; 2 w'm 


e 



d it)*. 


(9) 


where < Rj > and < 8*3 > are averages of Rj and B*j over the range of integration. This 
equation is reduced further, since for plane wave modes Rj and Bjj are independent of cj'. 
Taking these terms outside of the integral and integrating with respect to o)', the following 
limiting form of the intensity results for the case of low frequencies: 



M, (AP)* 


fl'* < R. > < B\ > 

Cj I |3 



-V4L'’ w'* 


w. 


( 10 ) 
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w. 
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V 


I 


1 


If, as is usually the case, the change in axial velocity and speed of sound across a stage 
is relatively small compared with the change in static pressure and swirl velocity, < R. > and 
<B,3 > can be much simplified. Using bars where appropriate to denote average values 
across the stage, the downstream values and are as follows: 






-O+M,) 

A? — 

2 + M. 

PpS* 


Since, for consideration of low frequency core noise, only the downstream tadiated 
wave is directly important, only the downstream intensity flux is written out in its approxi- 
mate form. Thus, 



(A P)^ 

2 a’/* 



(1 (1 +MJ* 



AP 



M. 
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L'_ 


1/ I ' * ** 

- Ljj (jj 
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— A P — 

or, since M* < < 1 and M, < < 1 for most turbine blade or vane rows, 

* - J * 

P« a* 
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Thf peak acoustic intensity can be determined by differentiating the integrand of equa* 
tion (9) with respect to oj 'and equating to zero. It is found to occur when cj ' -yfljL ’ . 
If this frequency satisfies the condition cj 'M^ < 1, the peak radiated acoustic intensity in d 
bandwidtli given by (tj ' — A cj cj + A w 1 is, so long as A co' < < cu' , 

ni m HI * 



( 12 ) 


P 



3 


From this approximate analysis, it can be inferred that for cases where the no'sc due to 
convccting temperature nuctuations is predominantly propagated in the plane wave mode, 
the acoustic intensity tlux within a fixed bandwidth is proportional to the turbine stage op- 
erating parameters (A P)^ - • ***'-' niean square intensity of the temperature fluctuations, 

and tlie transverse “correlation area” of these tluctuations. Tlie axial length scale does not, 
therefore, contribute in sucli cases to the peak intensity on a fixed bandwidth basis. Together 
with the axial velocity, however, it docs fix the frequency at which peak noise occurs. 

The range of validity of t* esc particular results can be demonstrated by calculating a 
range of cases using numerical integration of the trinic integral of equation (8). For exam- 
ple, the narrowband spectra of noise due to temperature tluctuations of fixed intensity and 
transverse length scale passing through the first row of exit nozzle guide vanes on a JT3D 
engine operating at a typical approach speed is shown in Figure 2 for a range of axial length 
scales. Also shown art some results for the low frequency approximation given by equation 
(12). It can be seen that the larger axial length scales generate approximately equal levels 
of peak noise at low frequencies and are reasonably well represented by the approximate form 
of equation (1 2). As the length scale decreases, noise is generated over a wider frequency 
range with slightly higlier peak values and the approximation is seen to be inaccurate. The 
changes in peak level and spectrum shape for axial length scales less than about 0.4 of a ra- 
dius are due to noise propagating in acoustic modes other than the plane wave mode. 

The effect of varying transverse length scale for constant intensity temperature fluctua- 
tions is shown in Figure ?. The peak one-third octave band power level is plotted aginst L' 
(which is taken equal to ^or comparison purposes) for three non-dimensional axial length 
scales at the iume operating conditions as used iit the previous figure. It can be seen that at 
the lower non-dimensional transverse length scales, the power level increases directly with 
the transverse correlation area and then deviates below this dependence at higher values of 
Ly and L\ Tliis effect is due to the propagation of higher circumferential modes and is seen 
to be slightly more pronounced for the smaller axial length scales which generate peak noise 
at higher frequencies. The dependence of peak one-third octave band power levels on the 
transverse correlation area of the temperature fluctuations is thus a reasonable approxima- 
tion over a reasonable range of transverse and axial length scales. As shown earlier, however, 
inaccuracies will be encountered if cases with transverse length scales on the order of a radius 
are considered and where in addition the axial length scale is less than 0.4 of a radius. 
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Figure 2 


Figure 3 
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spa'll ol' suuiui 
liair blailc chord 

coclTicicnt of cxpcnsion at constant volume 
coefficient of expansion at constant piessure 
number of rotor blades 

amplitude of acoustic pressure due to temperature fluctuations passing 

through turbine blade row 

drag force 

cyclic freiiuency 

lift force 

force vector 

non-dimensional spectrum of C (u', k'^. k'^) 

Heaviside function 

xTT 

acoustic intensity vector 

harmonic of inflow temperature distortion 

wave number vector . k^, . k, ) 

length scales (L^, of temperature fluctuations 

Mach number 

acoustic pressure 

steady state pressure 

representative radius of turbine stage 

cross-correlation function of temperature fluctuations 

radiation factor 

unsteady entropy 

time 

unsteady velocity vector (u, v, w) 
steady state velocity vector (U, V. W) 
rotor tip speed 
axial coordinate 
circumferential coordinate 
spanwisc coordinate 

Ui Ji 

ratio of specific heats 
delta function 
tagger mglc of blades 
density 
temperature 

frequency, radians per sec 
time displacement 
phase angle 
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Subscripts 

0 

I 

c 

i 

P 

X 

y 

z 


steady state 
upstream 
downstream 
convected wave 
integer 

acoustic (pressure) wave 
axial component 
circunU'erential component 
spanwise component 


Superscripts 


c 

P 

1 

2 

i 


convective wave 
pressure wave 
lower cutoff value 
upper cutoff value 
integer 

non-dimensional for dependent variables 
displacements of independent variables x', y', . 


Non-Dimensional Variables 



L/r 

I Jr 

^!0l 


car 
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SUPPLEMENT 

DERIVATION OF CONDITIONS FOR ^HE PROPAGATION 
OF ACOUSTIC. VORTICITY AND ENTROPY WAVES. AND THE 
JUMP CONDITIONS ACROSS /VN ACTUATOR DISC 

Using the notation illustrated in Figure I where small unsteady velocity components 
u = (u. V, w ) are superimposed on a steady and constant velocity field U = ( U. V, W ) relative 
to the stationary coordinate system x = (x. y. z) . the linearized equations for continuity, 
momentum (Navier-Stokes) and energy are. neglecting viscous and heat conduction effects: 


+ (U.V) p+p^ (V.u) = 0, 

(Al) 

n du 

(A2) 

+ (U.V)s = 0. 
dt 

(A3) 


In these equations, p. p and s are the perturbation pressure, density and entropy respectively, 
and and p^F are the ambient density and fluctuating force per unit volume appi ed at the 
cascade. In addition, the linearized eiiuations of state and second law of thermodv namics can 
be written as: 




= 0 , 


(A4) 


S + — p - -B-O =0, 
p 0 

O 


(A5) 


where P is the ambient pressure. 0^^ the ambient temperature R the gas c. nslant the coef- 
ficient of thermal expansi-'n at constant pressure and 0 the perturbation tempeiature. 

If the rotor is represented by an actuator disc at the x = 0 plane, and if the How \ariables 
arc assumed to have tlie form p (x. t) = p (k. cj) e’^- - etc., equations ( A1 ) - 1 A.>) become , 
away frorri the rotor plane x = 0, (where F = 0), and dividing by e'*-- 


(U . k - cj)p + Po (u . k) = 0. 


(U . k - w) u +_^ p = 0, 
Po 

(U . k - w) s = 0. 


(A6) 


Equations (A4).ind (A5) remain of the san-e form where, without ambiguity, all dependent 
variables arc functions in wave number/frequency space. 
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The seven equations given by (A4) to A6) represent u linear set of homogeneous equations 
in seven unknowns. The condition for a non-trivial solution of tliis set of equations is that the 
determinant formed by a matri.x representation of these equations be zero. This condition is 
satisfied when 


(U.k-w)^ Ik P - (U . = 0. 


(A7) 


where ( 2? ) is replaced by the speed of sound a. 
f>o 

Inspection of this equation shows it to be a 5th order polynomial in k^^ , one root of which 
is a triple root given by U . - cj = 0. or 


Vi 


(c) = cj - Vk^ - Wk 
S 2 L. 

U 


(A8) 


where the superscript (c) is used because it can be seen that the form given for k^^ gives rise to 
a wave convecting downstream with the mean flow. Substitution of this form for k^ into equa- 
tions ( A4 ) - ( .Ah) yields the tollowing relations lor variables convecting with the mainstream: 


Uck,<'> + v^k^ + w^k. 


e 

c 

Pc 


= 0 
= 0 
= s. 


0 Sc 

7C, 

PQ S 

c ' 

yc I 


(A9) 


where represents any constant entropy waves that may exist convecting downstream to the 
rotor plahe. 


Waves of this type represent vorticily and entropy waves propagating with no associated 
pressure fluctiKuions and w'ith phase velocities equal to the component of the mainstream vel- 
ocity in a direction normal to the wave front. The wave ilenneii by equations in ( A9) thus 
represents a spatial vorticity and entropy distribution that is convected downstream with the 
mainstream flow. 
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The remaining two roe’s give rise to upstcam and downstream acoustic waves (denoted 
by subscripts I and 2 resix'ctively). 


l(p) 

»i 


n, -d +k^) 

(1 -M[ ) 


V (P> 


->/n! -o-M.;)<kJ +k^) 


(1 -M,^) 

*2 


where k^. +M^k^ - cj/a. 


and I 


-(Uj.V..W.) = (M, ,M .M ) for i=l,2. 

t i 


Substitution of the results in equations (A4) - (A6) yield the following relationship that 
relate all acoustic perturbation variables in terms of acoustic pressure: 

_l((p) 


w. 


p, (y. .k -co) • 
_!l2___Pp 

~K 

* p 

p„(yi k-w) 

1 

IT Pp, 


(All) 


T (^) 


where i = 1 and 2 refer to upstream and downstream regions respectively. 
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It has been shown that in the absence of conductivity, the unsteady flow field away from 
the rotor cascade, can be represented by upstream and downstream propagating acoustic waves 
and vorticity and entropy waves conveeting with the mainstream How. In order to see how 
sound can be generated as \-.>i ticity and entropy waves interact with the rotor, it is necessary V 

to satisfy steady and unsteady jump conditions across the rotor disc, f luctuations in force 
are experienced at the rotor due to the unsteady flow causing variations in lift on the blades 
and this gives rise to one source of acoustic waves. Additional noise is generated, however 
due to entropy and vorticity waves experiencing a rapid change in static pressure across the 
rotor. The appropriate jump conditions are now derived and will be applied to the present 
problem where convected and acoustic waves are present both upstream and downstream of 
the rotor. 

Steady Flow Jump Conditions 

Using the same notation as in Figure 1. with the cascade plane being represented by an 
actuator disc at x = 0, and with upstream and downstream How variables represented by sub- 
scripts 1 and 2 respc(.tively. The steady equations of continuity and momentum for a How 
through the actuator disc are; 


V.(PoU) = 0, 


(A12) 


Po(MY)U + VP =p„F^6(x) (A13) 

where P^F^is now the steady force per unit area at the actuator disc. 

Integrating equation (A1 2) across the disc from - e to e and considering the limit as 
c -*• o, it is easily seen that only terms that contain variables differentiated by x need to be 
retained. Thus 


c 

lim r (p^U)dx = 0, 

€-»o J 3x 


where upon 

(At4) 


Similarly, the x component of equation (A1 3) becomes 



^jdx 

3xj 


Po COS? + sin t) 


(A15) 
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where and represent the steady lift and the drag at the cascade respectively. The lift 
is taken in the direction normal to the airfoils of the cascade and , the drag, is taken in the 
direction parallel to the airfoils. Also the bar notation as in ^ represents the average of the 
upstream and downstream values i e."^ = '/i (p^ * Po ^ 

In order to integrate across the singularity at the rotor plane, it is >.xpcdicnt to express 
both the upstream and downstream variables within one function by use of the Heaviside 
(or >itep) function. In this way, the axial mean stream velocity and its derivative are repre- 
sented throughout all space resivcivcly as 


U(x) = U, H(x)+ Uj |1 -IUx)l, 

= (U, -Uj) = AU 5(x), 

ox dx 


where (Uj - Uj ) is replaced by AU and the fact that ‘he differentiai of a Heaviside function 
with its argument is a delta function of that argument is applied. This technique is applied 
throughout this appendix, but in particular, equation (A15) becomes: 


Pj,U AU + AP = - p„ (F^ cos 5 + D„ sin ^). 


(A16) 


Similarly the y and z components of (A13) reduce to 


p U AV - -p (F sin ^ + D cos^l. 

(A17) 

AW = 0 

(A18) 


Unsteady Flow Jump Conditions 

Integrating the unsteady lineraized equation of continuity as given by equations (A1 ) 
across the x = 0 plane from - c to e yields, as before. 


lim 
€ -♦ o 



(pU + P„ u) dx 


0 . 


i.e. p,U, p,Uj +p„^u, 


= 0 . 
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The X component of the unsteady momentum equation, equation (A2) is integrated in 
the same way across the actuator disc. 


“it 




Si-/ 


p (F COS { • D sin () 6tx) dx 


+ / p(F„ cos^-D„ sinO^(x)Jx 


where p^F and p^D represents the fluctuating lift and drag per unit area at the cascade in the 
directions normal to and along the airfoil respectively (as shown in Figure 2). 

Using the same methods as with the steady flow jump conditions, this equation can be 
reorganized so that 

J |pUAU6(x)+ p^u A U8 (x)|.dx + {p„Uu)f^ - J | up^ AU + uUAp^|5(> 


+ IP)4 Po cos ^ - D sin {) + p (F_^ cos ? ~ sin ^), 


whereupon. 


pUAU + p^ ^ U, u, - p^ ^ Uj Uj - uU A p„ + Ap 


= p^ (F cos 5 - D sin {) + P (Fq ^ 


Tlic y and z components of the unsteady momentum equations can be considered in a similar 
manner reducing with help from the steady jump conditions to: 

AV (pU + P_ u) + p^ U, V, - p^ Uj Vj - p^ AU - vU Ap„ - P„ (F sin f + 0 cos 

‘ ' U 

- p (F^ sin ? + cos ^), 

Po|U, Wj -Po^UjWj -p^AU - Uw Ap^ =0. (/ 


'V 





Application of Jump Conditions Across the Rotor Disc 

It is assumed that vorticity and entropy waves convecting with the mainstream exist both 
upstream and downstream of tlie rotor and that the only source of acoustic waves is at the 
rotor. If the subscripts p and c for the unsteady variables indicate tliat the variables are due 
to acoustic pressure waves or convected waves respectively, and if as before the subscripts I 
and 2 denote upstream and downstream respectively, the upstream and downstream variables 
take on the form 


P, C, 

p. = p + p 

"J ^Pj 

with similar expressions for the other variables. 

Substitution of these expressions in the unsteady jump condition for continuity (A19); 

+ + u^^)-p„^(Up^+ u^^) = O. 


Applying the relations in (A9) and A1 1) this equation can be written relating the upstream 
and downstream pressure waves as follows; 



where = U.k^''’ - cu 

Further reduction using the steady continuity conditions and reorganizing, 


I 



k‘^n 

X . 


fu. 

Cl 




1 

<X 

Pp.- 

J 

_ ^ 

p = P u _ P U 

(A23) 


In like manner, the jump condition for the x-component of momentum, equation (A20) 
becomes, using the notation introduced earlier in this appendix: 
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/2( U,(Pp_ + p^^) + Uj(Pp^ + p^^)l AU + p„^U,(Up^ + u^^)-p„^Uj(Up^ + u^,^) 

'/2lU,(«p^ + u,^) + U,(u^^+u^.^)lAPo 

+ A Pp = ^ (F cos ^ - D sin 0 + ("^ +^) (F^ cos ^ sin ^) 


Using the equations (A9 and (A1 1) as before, together with the jump condition for tlie 
x-component of the steady momentum equation, it can be shown that this equation reduces 
to: 


>iAP >4Po(AU)’ KV/’! 




■/.p.cAu)’ p;u,k''” 




= Pp (F cos ^ - D sin {) - A P _L - ^U, u^ + p; Uj u^ 


(A24) 


The jump condition for the y-component of momentum (equation A23) can be obtained 
using similar methods. Thus 


'/2p„ AU kj < t 

./,AV{— ~-^\p U, + p U I 

(^af “i j Po} °> ' °2 ■/ '’ 


(■/2PoAU k, ) 
i/j A V.’ ^ li- 


2 Oj 


= p ( F sin ^ + D cos {) - '/i AV ( p u + p u ] - Vi j p„ U + p U ( (v - v ) 

(A25) 

and the z-component of momentum (equation A22) becomes: 


^o,“l ■■ ^ 0/2 


P = W - W 

« D C . C 


(A26) 



Now, the term u can be eliminated between equations (A23) and (A24) yielding an equation 
relating the up and downstream acoustic waves of the form 


Pp, +A,jP =p„(FcosS-Dsin{)-APjL _ ~ (u] (A27) 

I 2 o 1 2 

where 


A„ = I 


SiAP (AU)’ 


P a. 

'^o I 


J. 

^ a^ 

I 


Po U .«, 


(U, - U ) - 


P a* 
"2 ‘ 


( 'MP ^ 
^12 = - - 


(AU) 


P a^ 
2 


P aj 

2 


P aj 
°2 ' 


A second relationship between the up and downstream pressure waves can be obtained 

by substituting for u ,v and w from equations (A23) (A25) and (A2o) respectively m 
‘■’2 *^2 '■'2 

the second equation from (A9). Thus: 


Aj, Pp + Ajj Pp = p„ (F sin ^ + D cos ^ 


(c) 


M/2(Po^U,+P„^U,)^ APo +Po,^ AV] u^,_ 


(A28) 


where 


Vik 


(C) 


U, K':'i '-C, \ 


A:. = — ^ 

Poj ' hi ' ] ( ^ ’ 




- '/2(Po u +p U ) 


p„.«. 


u. 


'22 


(Po ^2 +Po U,> { - 

D 1 2 J r,2 


k <>'’ 


a- 


'•^AVk^ P,,^ U, 
To ^2 




^ ^ (k^ + K;). 

Poj«2 
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So long as the deteuninant (A,, Ajj - A^^ A^, ) is non-zero, a solution for the upstream 
and downstream pressure waves can be obtained in the form: 




B 


4i 




.i= 1.2 


(A29) 


where the non-dimensional coefficients are readily determined and the normalized up- 
stream temperature fluctuation convecting with the main stream has replaced s^ by using 

equation (A9). 7C^ 


It can be seen that acoustic waves are generated not only by fluctuating lift and drag at 
the rotor, but also by entropy (temperature) waves and the axial component of inflow velocity 
distortion as they experience a change in mean streani flow conditions across the rotor. The 
fluctuating lift and drag forces are moving with the blades of the rotor and generally generate 
high and low frequency noise depending on the velocity fluctuations giving rise to tlie nuctua- 
ting angle of attack at the blades. The noise generated by temperature and axial velocity 
fluctuations, which is the main concern of this paper, will generally be of low frequency and 
only significant when large changes in pressure, density and mean axial velocity are ex- 
perienced across the rotor. 

Upstream and Downstream Radiated Acoustic Intensity Flux 

Once the Fourier transform of a pressure wave p (w, ,kj ) is known in a constant 
moving medium, the acoustic intensity flux vector can be derived by applying the work of 
Bloklintsev. Thus, 


|p(co, k , kjl^ 

I_(w,k , k )= 

Po a 


(U + an) 
a - y . 


where J[ is a vector giving the magnitude and direction of acoustic intensity at frequency to 
and in wave numbers k and k , and where n is the vector normal to the wave front, and 

y z u , 

where the pressure waves have space-time dependence of the form e ' c ). 

For a typical rotor situation where the mean radial velocity is essentially zero (except of 
course for leaned blades or vanes) the upstream and downstream intensity flux per unit area 
radiated by a pressure wave is calculated from 


|p(to,k ,k,)| 

I, (co,k,,k^) = R,(to,k.k) i = 1. 
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where the non-dimensional R. (u.k^ , k^) is ^ven by: 


Ri=(-1) 


I 


(P) 


+ 


* r~^ — — r 

yJK. + k; + k 


1- 


CJ 

a, 


(A30) 


In order to derive the acoustic inteiisity llux at a given freciuency a. id for al! wave numbers, 
I and I have to be summed over all aj^plicable k and k and the solution doubled to 

*2 y ^ 

account for an additi(,nal set of pressure waves of the form p = p (-k , -k , -co) e ‘ * 


The applicable range of k^ and k^ correspond to those wave numbers that give rise to propaga- 
ting as opposed to decaying acoustic waves. The condition for propagation is that 


k^^^’and k,^'’’ are real. Thus, the upstream and downstream acoustic intensity flux per unit 
area at frequency cj is calculated from: 



where 




j 



-I + ( 1 - M. 


^ > (v) 


( I - M ■ 


) 


(A31) 


k<« = 


V 


(-ly 

r 


(^) 


1 - M - M 


j.i = 1,2. 
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The numerical calculation of these integrals is reasonably straightforward except for the 
fact that the determinant (A, , A,, - A^ A,j ) goes to zero at the limits of integration. In 
this event, the upstreair and downstream radiating pressure tend to infinity giving nse to a 
singularity in the integral of the form 1/x which is non-integrable. This limitation in the 
analysis is common to many solutions of linear systems that do not contain damping due to 
either dissipative or non-linear terms. Recognizing the fact that these infinities do not occur 
in reality , and that in many cases of interest, the bulk of the radiated energy is concentrated 
in modes not immediately at or next to cutoff, no serious error is introduced in the inte- 
gration. 
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APPENDIX B-2 

THERMOCOUPLE VOLTAGE SPECTRA 

One-third octave band spectrum levels of compensated thermo- 
couple voltage are presented for thermocouple elements identi- 
fied by the following code: 


THERMOCOUPLE IDENTIFICATION CODE 



; 
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THIRD OCTAVE BAWD TC VOLT AOiUVlt-dB - i H I i ' THIRD OCTAVE BAND TC VOLTAOt LEVEL - rt 



































THIBD OCTAVE BANO TC VOLTAGE LEVEL^ 

































APPENDIX B-3 


THERMOCOUPLE CORRELOGRAMS 

Normalized auto and cross correlation function machine plots 
are presented in this Appendix. The relevant transducer loca- 
tions are illustrated on each plot. Selected cross-correlations 
between thermocouple elements and wall pressure transducers 
are also included. 


THERMOCOUPLE IDENTIFICATION CODE 
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